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ABSTRACT 



This thesis describes the development of computer algorithms for experimental 
measurement control and subsequent signal processing of transient signatures to syn- 
thesize scattering impulse responses of scale model targets. The theories behind tran- 
sient scattering are considered in order to construct the algorithms. Up-to-date 
hardware and software technology are selectively implemented to optimize the resultant 
signal-to-noise ratio. The detailed explanation of the hardware and software operation 
is provided. A noise model for the system is also discussed. Measurement and pro- 
gramming validations are described, where comparisons are made with numerical com- 
putations for transient scattering by selected canonical targets. Results are documented 
for futher research in the transient scattering problem. 



TABLE OF CONTENTS 



I. INTRODUCTION 1 

A. BACKGROUND 1 

B. OVERVIEW OF THIS EFFORT 2 

II. SYSTEM FACILITIES AND HARDWARE 4 

A. LABORATORY DESCRIPTION 4 

B. ANECHOIC CHAMBER 6 

C. HP54120T DIGITAL PROCESSING OSCILLOSCOPE 6 

D. GURU II HARDWARE 8 

III. THEORY OF TRANSIENT SCATTERING MEASUREMENTS 11 

A. FREE-FIELD RANGE SYSTEM REPRESENTATION 11 

B. MATHEMATICAL MODEL 13 

C. TRANSIENT RESPONSE SOLUTION 14 

IV. ACQUISITION AND SIGNAL PROCESSING ALGORITHMS 17 

A. GURU II INTERFACE 17 

1. GPIB Principal 17 

2. Basic Software Requirements 18 

B. ACQUISITION ALGORITHM 18 

C. DECONVOLUTION ALGORITHM 20 

D. EXAMPLE OUTPUT 25 

V. EXPERIMENTAL RESULTS 27 

A. CANONICAL MODEL SCATTERING 27 

B. SCALE MODEL SCATTERING 32 

C. SYSTEM NOISE ESTIMATION 33 

VI. CONCLUSIONS 43 

A. SUMMARY 43 

B. FUTURE CONSIDERATIONS 44 

iv 



APPENDIX A. THE IMPULSE RESPONSE DECONVOLUTION 



46 



APPENDIX B. THE ACQUISITION QUICK BASIC PROGRAM 67 

APPENDIX C. CHARACTERISTICS OF LOW LOSS CABLE 78 

LIST OF REFERENCES 79 

INITIAL DISTRIBUTION LIST 81 



v 



LIST OF TABLES 



Table 1. IIP54I20T CHANNEL CHARACTERISTICS 7 

Tabic 2. HP54120TTDR AND TIMEBASE CHARACTERISTICS 7 

Table 3. TYPICAL NOISE 8 

Table 4. DATA OUTPUT TRANSFER RATE AND DATA RECORD 

LENGTHS 8 

Table 5. FREQUENCY DOMAIN FUNCTIONS 11 

Tabic 6. MATHEMATICAL NOTATION 14 

Table 7. CALIBRATION SPHERE PARAMETERS 23 

Table 8. 90 DEGREE THICK WIRE PARAMETERS 26 

Tabic 9. 12.2 CM SPHERE PARAMETERS 28 

Table 10. 45 DEGREE THIN WIRE PARAMETERS 29 

Table II. 90 DEGREE THIN WIRE PARAMETERS 30 

Table 12. S.I CM SPHERE PARAMETERS 31 

Table 13. 15.5 CM SPHERE PARAMETERS 33 

Table 14. DIMENSIONS OF SCALE AIRCRAFT MODEL I 34 

Table 15. DIMENSIONS OF SCALE AIRCRAFT MODEL 2 35 

Table 16. DIMENSIONS OF SCALE AIRCRAFT MODEL 3 36 

Table 17. DIMENSIONS OF SCALE AIRCRAFT MODEL 4 37 

Table 18. STATISTICS OF NOISE WAVEFORM 39 

Table 19. STATISTICS OF NOISE MODEL 40 

Table 20. STA TISTICS OF THICK WIRE SIGNAL PLUS NOISE WAVEFORM 41 

Table 21. CHARACTERISTICS OF PORT RETURN CABLES 78 

Table 22. INSERTION LOSS 78 



VI 



LIST OF FIGURES 



Figure 1. General layout of the scattering range 5 

Figure 2. System representation of the TESL 12 

Figure 3. Deconvolution input waveform 21 

Figure 4. Overlay 23 

Figure 5. Subtraction waveform 24 

Figure 6. Thick wire measurement 25 

Figure 7. The 12.2 cm diameter sphere validation 28 

Figure 8. The 10 cm thin wire validation at 45 degrees 29 

Figure 9. The 10 cm thin wire validation at 90 degrees 30 

Figure 10. The 8.1 cm diameter sphere validation 31 

Figure 11. The 15.5 cm diameter sphere validation 32 

Figure 12. Vertical and horizontal wing scale aircraft model 1 34 

Figure 13. Scale aircraft model 1 at 000, 030, 090, 180 degree aspect angle 35 

Figure 14. Scale aircraft model comparisons at 000 degree aspect angle 36 

Figure 15. Scale aircraft model comparisons 090 degree aspect angle 37 

Figure 16. Noise waveform statistics 39 

Figure 17. Noise power 40 

Figure 18. Signal plus noise waveform statistics 41 



vii 



ACKNOWLEDGEMENTS 



I wish to thank Professor Michael A. Morgan for his thoughtful and patient guid- 
ance throughout the course of this research. 

This thesis is dedicated to my child and my wife, Oy, who for over two and a half 
years has kept the home fires burning while I have pursued my education. I could not 
have done it without her. 



vrn 



I. INTRODUCTION 



A. BACKGROUND 

In circuit theory, it is well known that the characteristics of a linear time-invariant 
system can be obtained by exciting the circuit with an impulse function. The theory of 
transient scattering measurements is analogous to this method, and it also provides a 
unique parametric description of the scatterer via the complex pole locations in the 
Laplace domain transfer function. There are two approaches to obtaining the impulse 
response. One utilizes frequency domain methods and the other employs time domain 
techniques. The former has to sweep frequencies over the entire band but the latter uses 
direct pulse excitation. As a result, the latter can save time and expense. Therefore, 
direct transient scattering measurements are a viable alternative to the more conven- 
tional CW method for obtaining impulse responses of targets. 

In 1965 Kennaugh and MofTatt [Ref. 1] investigated the properties of the impulse 
response waveforms backscattered from metallic targets. The approximations for tran- 
sient scattering and impulse responses they made formed the basis for "ramp-response" 
imaging. 

In 1967 DeLorenzo [Ref. 2] presented a system that could be used to obtain an ap- 
proximate measurement of the impulse response of scattering bodies. His major devel- 
opment was an indoor laboratory scattering range which facilitated the time domain 
analysis. He showed the agreement between his calculations and the measurement of the 
impulse response of a rectangular plate. 

Also in 1967, Murchison and Falk presented a hybrid system for measurement of 
impulse responses [Ref. 3j. This measurement system was separated into two frequency 
ranges. For frequencies below 4 Ghz, they used an impulse generator while for the fre- 
quencies above 4 Ghz, they used a swept frequency CW source. 

In 1971 Baum [Ref. 4] demonstrated theoretically that characteristic resonances can 
be obtained for scattering objects. His technique was called the Singularity Expansion 
Method (SEM). It was postulated that the resonances could be obtained from the 
transient scattering measurement for a complex object. 

In 1973 Tesche [Ref. 5] applied the SEM to to numerically evaluate complex reso- 
nances of simple metallic scattering objects. 
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In 1975 Blaricum and Mittra [Ref. 6] proposed a technique for extracting the poles 
and residues of a system directly from its transient response. Their technique provided 
a number of numerical advantages. 

In 1980 the original Transient Electromagnetic Scattering Laboratory (TESL) be- 
came operational at the Naval Postgraduate School (NPS). Hammond [Ref. 7] devel- 
oped this laboratory for inverse scattering research based upon synthesized ramp 
responses. This laboratory utilized an outside ground-plane range. This range was 
limited to symmetric scattering objects which could be "mirror imaged". 

In 1983 the free-field TESL began operation at NPS [Ref. 8]. This laboratory is lo- 
cated in Spanagel 535. The TESL facilitated research efforts in high quality transient 
scattering measurements [Ref. 9] used for radar target identification and characterization 
using complex pole natural resonances [Ref. 10]. 

Prior to the current undertaking, McDaniel [Ref. 11] showed good agreement be- 
tween the computation and the measurement of the impulse responses for a sphere and 
a thin wire. These measurements and signal processing were performed using the 
Tektronix hardware set up, as originally employed in the 1980 TESL. In 1987 a rev- 
olutionary new Hewlett Packard digital processing oscilloscope became available. In 
addition, computational hardware had evolved dramatically since the early 1980's. This 
thesis describes the implementation of this new technology to enhance the performance 
of transient scattering and signal processing at the NPS TESL. 

B. OVERVIEW OF THIS EFFORT 

The main objective of this thesis was to incorporate new measurement and compu- 
tational hardware and to design and implement new software to automate high-quality 
transient scattering measurements. In this thesis, much work has been done in the de- 
velopment of the transient scattering work station which has the ability to measure and 
synthesize the scattering impulse response of scale model targets. The major hardware 
includes a microcomputer, a digital programable oscilloscope (DPO) and a shielded 
anechoic chamber. Automatic transient scattering measurements have been imple- 
mented by utilizing an IEEE-bus controller. 

With the new HP54120T DPO, the impulse responses of scale models are synthe- 
sized at a high signal-to-noise ratio (SNR). The acquisition and signal processing algo- 
rithms run much faster on the IBM PC-AT vis-a-vis the Tektronix interpreted BASIC 
computer previously used. The programs now use QUICK BASIC and FORTRAN 
compilers. Validations of software were made by way of comparing experimental 
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measurements with computations for simple canonical targets. Various scale models 
were measured and documented in detail. Furthermore, a noise calculation was provided 
for further research. 

This thesis is broken into six chapters. Chapter II describes system facilities and 
hardware. This includes laboratory description, anechoic chamber, HP54120T digital 
processing oscilloscope and the GURU II hardware interface. 

Chapter III presents the theory of transient scattering measurements. This is com- 
posed of the free field range system representation, mathematical model and transient 
response solution. 

Chapter IV deals with acquisition algorithms and signal processing. This chapter 
covers operational description, GURU II software interface and the deconvolution al- 
gorithm. In addition, some selected outputs are illustrated. 

Chapter V emphasizes the experimental results. Two types of targets are presented. 
One is a canonical target and the other is a scale model target. The noise estimation of 
the system is also discussed in this chapter. 

Chapter VI summarizes all work and provides some future considerations. 
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II. SYSTEM FACILITIES AND HARDWARE 



A. LABORATORY DESCRIPTION 

Figure 1 on page 5 shows the general layout of the scattering range. It is composed 
of 

• Digital Programmable Oscilloscope (DPO) 

• Step pulse generator 

• Low loss cables 

• Broad band solid state power amplifier 

• Broad band transmitting and receiving antennas 

• Anechoic chamber 

• Personal computer (PC). 

Note that the DPO includes the sampling and pulse generator unit. The DPO is 
used as both the master pulse source and as the receiver. The antenna connection cables 
were obtained from the DPO accessory’ kit for use in low loss network measurements. 
They are flexible for hard-to-access areas. The characteristics are shown in Table 21 and 
Table 22 on page 78 in Appendix C. 

Transmitting and receiving antennas are of the double-ridged horn type. Their fre- 
quency range is 1-12.4 GHz, with less than IdB of ripple in the boresight gain over this 
multi-octive band. The Avantek solid state power amplifier has an advertised 3dB 
bandwidth of 2 to 6 GHz but, by our observations, has a usable bandwidth of 1 to 7 
GHz. A 15 V supply is needed for the amplifier. 

The anechoic chamber serves two purposes. First, it is used to isolate the target 
signal from the electromagnetic noise and interference outside of the chamber. Secondly, 
the absorbing material on the walls, floor and ceiling of the chamber are designed to 
provide minimal reflections, thus approximating unbounded space. The computer, 
which serves both acquisition and signal processing rolls, is an IBM-AT with 640 K 
RAM, a 20 megabyte hard disk, and a GPIB-PC2A interface card. The GPIB-PC2A 
was added to the PC in order to have the ability to talk or listen to the HP54120T DPO. 
The PC was set to run at 9 MHz. It has two floppy disk drives, and two output ports. 
One is a parallel port for the printer and the other is a serial port for the plotter. The 
color graphic board is NEC GB1. The display is NEC EGA color monitor. 
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Figure 1. General layout of the scattering range 
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B. ANECHOIC CHAMBER 



The anechoic chamber is located in room 535 of Spanagel Hall. This chamber in- 
corporates a metallically shielded anechoic enclosure to minimize outside EM interfer- 
ence. Chamber dimensions are 20 ft. long with a 10 ft. width and 10 ft. height. Because 
targets are suspended within the chamber, there is no requirement that they have mirror 
symmetry, as is the case for ground plane based ranges. In the measurements that were 
taken for this thesis, the target is mounted on a thin styrofoam column at a distance of 
8.5 feet from the transmitting and receiving antennas. Both receiving and transmitting 
antennas are placed close together to closely approximate monostatic measurements. 
The chamber interior is covered by a special carbon impregnated foam absorbing mate- 
rial. The material is cut into longitudinal wedges along the side-walls, floor and ceiling 
to act as conduits of EM energy towards the back wall. The back wall is covered by 18 
inch long pyramids to act as slow absorbing transitions for incident fields. Further in- 
formation about this free-field TESL can be found from Refs. 9, 11 and 12 

C. HP54120T DIGITAL PROCESSING OSCILLOSCOPE 

The Hewlett Packard 54120T is a fully programmable 20 GHz, four channel, 
digitizing oscilloscope with a nine inch color display. It is capable of automated meas- 
urements, digital storage, and TDR measurements. This DPO has only been in pro- 
duction since mid- 1987 and represents the state-of-the-art in sampling technology. 
Previous to this thesis effort, a 6 year old 12 GHz DPO from Tektronix was used to 
perform transient scattering measurements in the TESL. 

The HP-DPO is comprised of two major components. These are the HP54120A 
Digitizing Oscilloscope Mainframe and the HP 54121A Four Channel Test Set. It also 
contains extensive self-tests to ensure proper functioning. Both HP 54120A and HP 
54121A were calibrated together as a system called the HP 54120T. The DPO can be 
remotely programmed by an HP-IB (IEEE-488) controller interface. The HP-IB com- 
plies with the IEEE 488.2 standard. The essential characteristics of the HP 54120T are 
summarized in Table 1 and Table 2 on page 7. 

It should be noted that the HP 54121A 4 Channel Test Set is sensitive to 
electrostatic discharge (ESD). Precautions should be taken to keep ESD from damaging 
the equipment. Details can be found in Ref. 13. 

Waveforms that are stored in the DPO memory can be sent to either an HP Thinkjet 
printer or an HP Plotter, both of which are accessed via the HP-IB bus. There are two 
ways to operate the instrument. One is manual operation via the front panel and the 
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Table 1. HP54120T CHANNEL CHARACTERISTICS 



Bandwidth 


20/12.4 GHz 


Number of channels 


4 


Full scale 


8 division 


Maximum scale 


80 mV/division 


Minimum scale 


1 mV/division 


Programable dc offset 


± 500 mV 


Input maximum 


±2V dc + ac peak 


Norminal impedance 


50 ohm 


Connectors 


3.5 mm 


Percent reflection 


less than 5% for 30 ps risetime 



Table 2. HP54120T TDR AND TIMEBASE CHARACTERISTICS 



Risetime 


less than 45 psec at 12.4 GHz bandwidth in average 
mode 


Timebase full-scale 


10 division 


Minimum scale factor 


10 psec/div 


Maximum scale factor 


1 sec/div 


Time interval accuracy 


less than 10 psec + 0.1 % of reading 



other is to operate by sending HP-IB bus commands via the computer. A major effort 
in this thesis was developing programs to provide easy operator control and signal ac- 
quisition of the DPO using the computer. These algorithms are explained in detail in 
chapter IV. 

Ensemble waveform averaging has been implemented in acquiring scattering signa- 
tures. The averaging reduces noise by 1 lJ~N , where N is the number of averages. The 
minimum discernible signal level in the digital processing oscilloscope system is approx- 
imately 35/t V. Typical noise amplitudes for the DPO sampling unit are shown in Table 
3 on page 8. 

Data display resolution is 1024 points horizontally x 256 points vertically (for dis- 
play times of at least 200 ps/div). The display can be selected in persistence or average 
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Table 3. TYPICAL NOISE 



20 GHz bandwidth, Avg = 1 


Noise 1.2 mV (RMS) 


20 GHz bandwidth, Avg = 256 


Noise 80 m V (RMS) 


12.4 GHz bandwidth, Avg = 1 


Noise 500 m V (RMS) 


12.4 GHz bandwidth, Avg = 256 


Noise 35m V (RMS) 


12.4 GHz bandwidth, persistence 


Noise 400 m F( RMS) 



mode and the number of averages can be specified as a power of 2, up to 2048. The data 
output transfer rate and the data record lengths are shown in the following table. 



Table 4. DATA OUTPUT TRANSFER RATE AND DATA RECORD LENGTHS 



10 psec/div < time/div < 20 psec/div 


100, or 400 points/record 


20 psec/div < time/div < 50 psec/div 


100, 400 or 800 points/record 


50 psec/div < time/div < 200 psec/div 


100, 500 or 1000 points'record 


200 psec/div < time/div < 1 sec/div 


128, 256, 500, 512 or 1024 points/record 


Data output transfer rate 


115 kbytes/s 



For the time domain reflectometer (TDR) system, the risetime is an adjustable value 
based upon the timebase setting. The minimum is 10 ps or .08x time/div, whichever is 
greater. The maximum is 5 x time/div. Usually channel 1 is used for reflection meas- 
urements, however, in this thesis, it was used as the source generator. 

The digitizer uses a 12-bit successive approximation AD converter. The full-scale 
range of the A/D is 640 mV. The least significant bit (LSB) of the A/D converter equals 
250 /xV. This gives one part in 2560, or more than 11 bits of resolution. Averaging can 
stretch out the resolution to 32 ixV. This improves the resolution to around 14 bits. 

The signal is sampled and digitized at a rate determined by the trigger rate, repe- 
tition rate, time base range, display mode, and number of channels turned on. Further 
information concerning the characteristics of the HP 54120T can be found from Ref. 
14. 



D. GURU II HARDWARE 

The Textronix supplied General-purpose interface bus User's Resource Utility 
(GURU II) is a combination hardware and software package that allows the use of any 
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of the IBM personal computers (PC), and some PC compatibles, to control GPIB 
programable instruments. In this section, the focus is on the hardware. The software 
is explored in the next chapter. 

The hardware consists of the PC2A GPIB interface card that plugs into the PC main 
board, and the GPIB cable. The GPIB PC2A interface card is an 8-bit half-size board 
which performs as a talker, listener and controller. It also can be implemented with the 
full range of talker, listener, serial poll, service request, and remote programming func- 
tions for as many as 14 devices on the bus. It can perform as a complete controller. 
Using the PC2A, data transfer between the system memory and the GPIB are performed 
at rates of more than 300 kilobytes per second. The PC2A is compatible with all revision 
levels of the IEEE-488 standard including the HP-IB. 

The GPIB-PC2A interface card consists of these major sections: 

• address decoding 

• buffering and data routing 

• interrupt arbitration 

• direct memory access (DMA) arbitation. 

• configuration switches and jumpers 

• GPIB-adapter-TLC (talker/listener/controller) 

• time of day clock with battery back-up. 

The address decoding monitors the address lines to recognize when the GPIB-PC 
I/O address is present on the computer I/O channel and enables read and write access 
to the GPIB adapter register. The buffering and data routing handles data transfer be- 
tween the IBM PC-AT I/O channel and the GPIB adapter through a bidirectional 
internal data bus. 

Interrupt arbitation recognizes when interrupts have been enabled or disabled and 
passes or inhibits them accordingly. DMA arbitation recognizes when DMA operations 
are enabled or disabled, and when the last transfer has taken place. It also routes the 
DMA request and acknowledges signals to the selected DMA channel. The DMA is 
used to transfer data directly between a peripheral and memory without CPU. 

The IEEE-bus controller forms an interface between the computer and the 
IEEE-488 bus. Conceptually, the IEEE-488 bus behaves in a way very similar to the 
system bus. Hewlett-Packard originally devised it [Ref. 15] and intended it to link to- 
gether programmable instruments in a laboratory or industrial environment. By con- 
trolling test equipment and measuring devices from the IEEE-488 bus, implementing an 
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automatic testing station is possible. A system under test is connected to the test 
equipment and measuring devices. The computer configures all equipment via the 
IEEE-488 bus and then reads the test results from the same bus. The GPIB has the 
following limitations: 

• Half of the devices attached to the GPIB must have power on. 

• The total length of the bus cable must not exceed 20 Meters. 

• No more than 15 devices may be attached to one bus. 

More detailed information can be found in Ref. 16. 
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III. THEORY OF TRANSIENT SCATTERING MEASUREMENTS 



A. FREE-FIELD RANGE SYSTEM REPRESENTATION 

The physical model of the system free field scattering range was first presented by 
Morgan in Ref. 12. The system block diagram is shown in Figure 2 on page 12. Im- 
portant components and parameters of the system are summarized in the frequency do- 
main representation in the following table. 



Table 5. FREQUENCY DOMAIN FUNCTIONS 



Pulse generator, driving the transmitting antenna 


*00 


Transmitting antenna response 


HAf) 


Receiving antenna response 


h r (J) 


The multiple scattering between target and chamber wall 


Hs(J) - - H c 


Direct couple of receiving and transmitting antenna 


HJf) 


The multiple scattering between calibration sphere and 
absorber 


H c scif) 


The multiple scattering between target and absorber 


mdj) 


The calibration sphere transfer function 


Hs(f) 


The target transfer function 


HKf) 


Thermal noise 


m 



The transient pulse source is the step generator in the TDR system. This generator 
not only supplies the low-level transmitted pulse, before amplification, but also acts to 
coherently trigger the sampling circuits of the DPO receiver. It thus is akin to the 
"master oscillator" in a coherent pulsed CW radar system. The risetime of the step 
generator is approximately 25 psec. This step pulse is first amplified and filtered by a 
broad-band GaAs FET amplifier. 

The amplified pulse, which no longer resembles a step waveform, is then fed through 
the transmitting antenna which radiates into the anechoic chamber. The radiated field 
emanating from the antenna is approximately an angular weighted spherical wave. At 
the target position the incident field is a quasi-plane wave. Because of the close prox- 
imity of the transmitting and receiving antennas, some strong coupling of energy is in- 
duced in the receiving antenna. 
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Figure 2. System representation of the TESL 
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This undesirable signal can be seen clearly if the time delay setting on the DPO is ad- 
justed to begin many nanoseconds before the target response time-window. Therefore, 
the time setting delay must be checked when the target is present. 

The initial calibration of the beginning of the target time-window was performed by 
using a rectangular copper plate, 1 ft wide and 1.5 ft high, as the test scatterer. The back 
scattering was observed with the copper plate turned broadside to the centroid of the 
horn positions and placed where the target would be mounted. The copper plate scat- 
tering signature was so large that it could be clearly discerned. The results of this 
measurement indicated a DPO time delay setting of 38 nsec was needed. 

Multiple scattering between the target and the chamber absorber is represented by 
the two way interaction arrow connecting HJiJ) to H s {f). The total noise of the system 
N(J), is due to both thermal emission from the chamber and receiver noise in the sam- 
pling front end of the DPO. These sources are combined together and denoted as N n (J) 
for the n-tli measurement that is made. To understand the theory of transient scattering 
measurements adequately, the mathematical model will be presented in the next section. 

B. MATHEMATICAL MODEL 

The five key parameters in Table 6 on page 14 play a major role in both the physical 
meaning and mathematical sense of the acquisition and signal processing algorithms. 
These five parameters can be related to the system representation, described in the pre- 
vious chapter, by the following three equations: 

r,W = HjV) //*(/) m [IIaV) + Uc<J) + Us(f) + U C sc W] + AT,(/) (3.1) 

nw = n,(f) Hr<J) m W A d) + IlcWl + a \(f) (3-2) 

1 jW = Hi(f) UrV) m ["a(D + 2/cW + Hl<J) + ulcV)] + A'jW- (3-3) 

These three equations represent the three respective measurements that are neces- 
sary for processing the scattering signature of a given target. As shown below, these 
measurements are for a calibration sphere, whose measured scattering can be compared 
to an accurate calculation; the "background" which denotes no target present, and the 
actual target scattering measurement. 

Equation 3.1 states that the response to the amplified pulse generator, X(f), due to 
the calibration sphere and the chamber is equal to the product of the pulse excitation, 
transmitting antenna transfer function //,(/), receiving antenna transfer junction //«(/) 



13 



Table 6. MATHEMATICAL NOTATION 



The calibration measurement waveform 


>'.(/) 


The background measurement waveform 


>':(/) 


The target measurement waveform 


Y 3 (J) 


The subtracted calibration waveform 


YJJ) 


The subtracted target waveform 


Y S V) 



times the sum of the antenna coupling HJJ) plus the calibration sphere 7/f(/) and the 
multipath multiple scattering Added to this product is the thermal noise NJJ). 

Equation 3.2 represents the background measurement, with no target present. The 
system response is similar to that previously considered, but with the absence of the 
transfer function for the calibration sphere. Equation 3.3 expresses the target measure- 
ment in the frequency domain. The intepretation follows the same pattern as equation 
3.1 except the target transfer function, /■/£(/), appears instead of that for the calibration 
sphere. The multipath multiple scattering in equation 3.3 also changes and is denoted 
by 7/ic(/). 

The three equations as stated above are three crucial measurements in the transient 
scattering problem. The backscatter measurement with the canonical (calibration) ob- 
ject }',(/), measurements without target (background) Y 2 (J) and measurements with the 
desired target Y } (J) are obtained directly using the time domain acquisition program in 
Appendix 13. This program was written in BASIC as part of this thesis effort. To speed 
up the execution of the program by a factor of about fifty times vis-a-vis direct inter- 
pretation, it has been compiled using the Microsoft Quick Basic Compiler. 

After acquiring accurate data, the next step is to process the data and extract the 
impulse response. This is the topic of the next section. 

C. TRANSIENT RESPONSE SOLUTION 

The three measurements that were just considered form the basis for estimating the 
"smoothed" impulse response of the target. The band-limited nature of the measure- 
ments does not allow synthesis of the true infinite bandwidth impulse response which is 
due to an incident delta function plane wave. Instead, the "smoothed" term means that 
we can estimate the transient target response due to some predefined incident signal (a 
double-Gaussian shaped waveform is used here) whose significant spectral content does 
not exceed the bandwidth of the measurements. Aside from the currently used solid- 
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state power amplifier, the bandpass range of the DPO-antenna system is about 1 GHz 
to 12 GHz. The spectral bottleneck is due to the GaAS FET amplifier, whose passband 
is about 1 GHz to 7 GHz. 

In order to synthesize the smoothed transfer function of the target //](/), the first 
step is to subtract the "clutter" portion of the responses in equations 3.1 and 3.3 which 
is embodied in equation 3.2. This subtraction results in equations 3.4 and 3.5 below: 



>M = >M ~ *M (3.4a) 

>M = Hjif) H r {J) XV) [H C S {J) + /&(/)] + AM (3.46) 

AM = ^(/) "AM (3.4c) 

>M = Y 3 V) ~ Y 2 {f) (3.5a) 

*M = HjV) HrV) m + //icOO] + AM (3.56) 

AM = AM - AM. (3.5c) 



Equations 3.4 and 3.5 are the very critical steps in obtaining the smoothed transient 
response of the target. The critical nature of the subtractions has to do with coherency 
in time of the subtracted clutter. The calibration sphere measurement }’,(/) had the 
background measurement V 2 (J ) subtracted in order to eliminate the clutter. The clutter 
signal is that which is present aside from the target and is composed of two main con- 
stituents. One is the direct coupling of the receiving and transmitting antennas, //<(/), 
and the other is direct scattering (as opposed to multipath) from the absorber in the 
chamber, // c (/). Since the measurements are made directly in the time domain, the sub- 
tractions are done in time domain. In a similar fashion, equation 3.5 can be obtained. 

The results, ) 4 (/) and }' 5 (/), are then transformed into frequency domain by discrete 
fast Fourier transform ( FFT ). This leads to the next step, which is to form an optimal 
deconvolution estimator for target response given by 



X 0 (f) "1(f) = 



m n't/) Ysdf) 

I >4(/)| 2 + C 



(3.6) 



where } 4 (/) is the complex conjugate of } 4 (/), Y sc {f) is the computed calibration sphere 
and C is the smoothing parameter. 
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(3.7) 



Yscd) = Us (f) 

x o(0 *- FFF -* x oU) (3.8) 

x Q {t) = C, c~ “> (f ~ '° )2 + C 2 e~ ^ (3.9) 

where, for a specified a, and a 2 in this "Double-Gaussian" ( DG ) pulse, the C, and C 2 
are chosen such that 

•*•<>(0=1 (3.10) 

and 



p+OO 

= 0. (3.11) 

Each of the two Gaussian waveform in equation 3.9 will be specified here by a "10% 
pulse width", defined by A/* such that 

(A t k ? 

e -*K 4 =0.1, (3.12a) 



which gives 



a k = 



4 

(A i k f 



ln(10) = 



9.21 

(A r k f 



(3.126) 



The estimator in equation 3.6 is due to Riad [Ref. 17] and results from an optimal 
least-square fidelity criterion for deconvolved target response. The result of equation 
3.6, X 0 iJ) HI , represents the frequency domain scattering response of the measured tar- 
get to an incident plane wave whose Fourier transform is X 0 (j). This incident waveform 
is specified by numerically computing the scattering response of the calibration sphere 
to it. This numerically generated response of the calibration sphere is embodied in 
Y sc (f) , as shown in equation 3.7. It is obtained from the SPRSCT program which uses 
the magnitude and phase data from the MIE series program. Both of these programs 
were written by Morgan. 

Both subtraction and optimal deconvolution algorithms are implemented in the 
Deconvolution Program in Appendix A. The procedure to operate this program will be 
explained in chapter IV. 
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IV. ACQUISITION AND SIGNAL PROCESSING ALGORITHMS 



A. GURU II INTERFACE 
1. GPIB Principal 

This section emphasizes the GURU II software interface. In order to under- 
stand the Acquisition program operation, the GPIB principle has to be understood. The 
General Purpose Interface Bus uses a twenty-four conductor cable. Sixteen of these are 
used to send the signal. Data is sent back and fourth between the listener and the talker 
in a bi-directional fashion. This means that both transmitted and received data travel 
on the same line. Eight lines are used for 8 bits of parallel data transfer. A byte is sent 
by serial mode but each bit is sent in parallel mode. Three lines act as the coodinator 
between the talker and the listener. These lines are assigned their duty as follows: 

• no data available line 

• data not ready line 

• data not accepted line. 

The talker on the bus controls the data available on the line and waits for the 
listener's response on the "not ready for data" line or "data not accepted" line. The re- 
maining five lines are the attention line, the interface clear line, the remote enable line, 
the service request line and the end or identify line. When the attention line is active, 
the talker and the listener will be specified. The controller has full access to five of these 
lines and has the capability to directly address other devices to command them to talk 
or listen. The interface clear line is used to set the turn-on condition. When the remote 
enable line is asserted, all devices on the bus are under computer control and nothing 
external can interfere. The service request line is used by a device to indicate a require- 
ment for attention and to request an interruption of the current sequence of events. 
There are two methods for the controller to determine the device and the message. One 
is the serial poll while the other is the parallel poll. The polls obtain numerical codes 
which show a specific event. The serial poll is a mechanism that the controller uses to 
obtain status from individual devices. As each device is polled, it returns a status byte 
which controller decodes to determine the device requesting service. The parallel poll is 
a mechanism for accepting and decoding the general status of as many as eight devices 
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simultaneously. The end or identify line is used by the talker to indicate the end of the 
message. 

2. Basic Software Requirements 

There are 3 files which the GURU requires in order to run the application pro- 
gram. These files are 

• GPIB.COM 

• IBCONF.EXE 

• CBIB.OBJ. 

The primary GPIB device drivers are read in via the CONFIG.SYS file at boot-up time. 
This file is in the root directory and must include the following DOS statement: 

• DEVICE= GPIB.COM. 

This causes the GPIB driver in GPIB.COM to be loaded when the system is booted. 

The IBCONF.EXE is the executable program to configure the GPIB.COM file 
containing the device drivers for each instrument. HPDEV is the name of the device 
driver for all Hewlet-Packard instruments, including the HP 54120T, whose decimal pri- 
mary address is 17. The IBCONF.EXE (Interface Bus Configure) program allows a 
match between the softwere and the hardware. The program prompts with menus of 
details about the hardware parameters associated with the interface card. This device 
for Hewlett-Packard is already configured in the system. When the IBCONF.EXE is 
run, the device drivers are mapped for each instrument in the display. By selecting the 
HPDEV, and editing it, the characteristics of the HPDEV show on the screen. This al- 
lows the operator to change the characteristics to be compatible with the hardware. 
After reconfiguration of the GPIB.COM file, the operator must return to the operating 
system and then re-boot the PC for the new parameters to be in effect. 

CBIB.OBJ is the object file which is linked with any compiled BASIC program 
which needs to interface with the GPIB.COM device drivers. In our case, the Acquisi- 
tion Program, written in GW-BASIC, is first compiled using the Microsoft QUICK 
BASIC compiler. This process creates an object (.OBJ) file. The CBIB.OBJ object file 
is then linked with the Acquisition object file. The final result is an executable (.EXE) 
file. 

B. ACQUISITION ALGORITHM 

The Acquisition program was written in Microsoft GW-BASIC. The hardware is 
interfaced with this program by way of the CBIB object file. This file comes with the 
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GURU II software. After compiling the source code, the object file of the Acquisition 
has to be linked to the CBIB object file in order to complete the creation of the Acqui- 
sition executable file. To run the program, the Brun3087 execution file is needed in the 
same DOS Directory. This file is supplied with Microsoft QUICK BASIC. The Ac- 
quisition program uses a friendly and interactive menu format. The program begins with 
a displayed welcome message, including the time and date. This screen also provides a 
brief procedure on how to measure the transient scattering response. The screen gives 
operator the choice of using the system default or the last response. The system default 
was designed for first time operators. The last response was designed for experienced 
operators. A main menu is displayed on the screen after the operator decides whether 
to use the system default or the last response. 

The main menu contains many importance parameters. A first item is the number 
of data points in the acquired time series. This number has to be a power of 2, but not 
more than 1024. This limitation could be modified in the program by defining a larger 
array size up to 2048 points. The second item is the number of sub-averages. This pa- 
rameter determines the number of values to average for each time point when the 
HP54120T is in the average mode. This parameter has to be an integer from 1 to 2048, 
also in a power of 2. The third item is the number of data blocks. This determines the 
number of sets of data in one record. The fourth and the fifth items are a user supplier 
identification for the target and the date of measurement. The sixth item tells the oper- 
ator about the typical sampling time interval. The seventh parameter is the type of the 
waveform being acquired on the measurement being made. As discussed in the previous 
chapter, there are three scattering waveforms which need to be acquired for each target 
being considered: calibration (sphere), background and actual target. The eighth item 
tells the operator the time window in nanoseconds. The ninth, tenth and eleventh items 
are the maximum vertical setting, the data file name output and the auto time setting. 
At the bottom of the menu, the operator is asked if he wants to change the items. 

The second menu appears on the screen when the operator decides to change some 
item on the main menu. This "overlay" menu contains all changeable parameters and 
allows a change to these parameters interactively. When the changes are completed, the 
operator can go back to the main menu by typing the letter "R". This brings the oper- 
ator to the main menu and displays the new item on the menu. When everything is 
ready, the operator selects "no change" to begin to automate the DPO. The program 
starts by initializing the DPO and setting all specified parameters as directed by the main 
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menu. The program stops when the DPO has finished the setting operation. The pro- 
gram then asks for the bandwidth of the measurement. There are two choices: 20 Ghz 
or 12.4 Ghz. To protect from errors in the initial setting, the program stops for the final 
check. As soon as the operator initiates the acqisition cycle, the beginning time shows 
on the screen. For a typical acquisition composed of 2048 averages and 1 data block, 
it takes about 7 minutes to finish the operation. In the meantime, the computer waits 
for the DPO to provide the raw data. When the DPO finishes this operation, the ending 
time is displayed and the PC opens a temporary hard-disk file for storing the raw data. 
Before the operation is ended, the header which is to be placed on the data file is dis- 
played on the screen. The operator can then observe any error in the data set, should 
there be any. 

The header is transferred as a separate file. The header provides important charac- 
teristics of the data, such as format type, number of points, number of averages, x in- 
crement, x origin, x reference, y increment, y origin, y reference and range setting. All 
of these parameters allow the program to recover meaningful data from the raw data. 
The remainder of the program constructs the standard data file structure and the hard 
copy output. 

At the end of the program, there is a pause with a query regarding the operator's 
intention to continue with an additional measurement. To complete the acquisition cy- 
cle for one target, the three measurements previously discussed have to be made in order. 
These three files will then be supplied to the deconvolution program in the next proc- 
essing step. For acquiring more than three measurements, the background must be 
measured for each target so it can be subtracted from the respective target measurement. 
The background measurement occurs either immediately before or after the target 
measurement. This reduces errors induced by long term nonstationary drifts in the sys- 
tem. The need for such temporal coherence of the target-background measurements has 
been verified experimentally during the course of this work. Typical outputs are as 
shown in Figure 3 on page 21. 

C. DECONVOLUTION ALGORITHM 

The deconvolution program was written in Ryan-McFarland (RM) FORTRAN in 
order to improve the run time. The source code is standard ANSI 77 FORTRAN so it 
may be recompiled under virtually any PC FORTRAN compiler such as Microsoft or 
Lahey. This program was originally written in Tektronix BASIC by Morgan [Ref 11] 
and was restricted to use on a Tektronix 4052 microcomputer. 
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DECONVOLUTION INPUT STEPS 



COMPUTATION SPHERE 



CALIBRATION SPHERE MEASUREMENT 




BACKGROUND MEASUREMENT 






Figure 3. Deconvolution input waveform 



The use of a standard FORTRAN-77 compiler generalizes the program so it can be 
run on any IBM PC compatible. This algorithm follows the mathematic deconvolution 
as described in chapter III. The program has graphic capabilities which are supplied by 
GRAFMATIC software. The procedure to compile the program is described in the RM 
FORTRAN manual. This program must be linked to the GRAFMATIC library in or- 
der to run properly. 

The program begins by displaying the standard input data file. This data File format 
has been standardized for use with all of the software support systems. The header of 
each data file contains the important information such as identification, type of meas- 
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uremcnt, date of measurement, time of experiment, the maximum time scale, data block 
and number of points. Types of measurement arc as explained in the previous section. 
The first data file needed for the deconvolution program is that for the computed scat- 
tering by the calibration sphere target. This computation data is obtained from the M 1 12 
scries program and SPRCT program, written by Morgan, as described in the previous 
chapter. The input data required for these programs are summarized in Table 7 on page 
23 . 

The sphere scattering computation is transformed into the frequency domain and 
displayed on the screen in order to observe the bandwidth. The bandwidth of the com- 
putation has to be approximately the same as the estimated bandwidth of the measure- 
ment. If the bandwidth is not appropriate, the operator should to terminate the 
program via a Ctl-C and run the MIE series followed by SPRCT program to provide 
appropriate data. 

In the next step, the program asks for the calibration and background file names. 
Both data are read in and plotted on the same axis in order to observe if time shifting 
has occurred, and if so, to be able to estimate how much. The operator is then allowed 
to shift one waveform either left or right and then re- observe the overlay ed waveforms. 
Once a satisfactory alignment is observed, the calibration waveform is subtracted from 
the background waveform and the difference waveform is transformed into the frequency 
domain. 

The program then stops and asks for the target file name. It then asks if the oper- 
ator desires to use the same background or another, depending on which background 
was closer to the target completion acquiring time. In a similar fashion, the target data 
and the background data are plotted on the same axis. If there is an observed drift, the 
program allows the operator to shift the target waveform to the left or to the right until 
both target and background arc lined up. Then the subtraction can be performed. The 
difference waveform is then transformed into the frequency domain. Typical scattering 
data being processed by deconvolution program is shown in Figure 3 on page 21. The 
overlay and shift is shown in Figure 4 on page 23 and the subtracted waveform is shown 
in Figure 5 on page 24. The target for this example was a 10 centimeter long "thick" 
copper wire, and the calibration target was a standard 0.1025 meter radius aluminum 
sphere. After the subtraction, the difference calibration waveform and the difference 
target waveform were tranformed into the frequency domain. 
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TARGET AND BACKGROUND MEASUREMENT COMPARISON 




Figure 4. Overlay 



Table 7. CALIBRATION SPHERE PARAMETERS 



Sphere radius in meters 


0.1025 


Distance from sphere in meters 


2.5 


Bistatic angle in degrees 


3.0 


Scattering plane 


E plane 


Time window in nanoseconds 


20.0 


Nyquist frequency in Ghz 


12.775 


No. frequency points 


257 
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SUBTRACTED WAVEFORM 




Figure 5. Subtraction waveform 

The three frequency domain waveforms were formed into Riad's optimal deconvo- 
lution estimator given in equation 3.6. The smoothing constant "C" is selected by the 
operator. This constant is set to optimize the denominator threshold level in the fre- 
quency domain. If it is set too large, oversmoothing will result which reduces the spec- 
tral content in both the low and high frequency regimes. If the threshold is set to zero, 
a "naive" deconvolution results. This provides vastly amplified noise in the signal and 
typically looks like "garbage". A typical value for C is 0.1. Riad's method yields an 
estimator for the scattered waveform of the target due to an incident plane wave signal 
having signal shape equal to that which illuminates the computed calibation sphere us- 
ing program SPRCT. This frequency domain estimator is transformed to the time do- 
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main and the result is displayed. At the end of the program, the operator may try' 
another smoothing parameter. If the result satisfies the operator, it can be stored in the 
file whose file name was provided by the operator. The program also allows the proc- 
essing of another target directly without rereading the computed data and the calibration 
data. 



D. EXAMPLE OUTPUT 

The output of the deconvolution program is in both the time domain and frequency 
domains, but the following example shows only the time domain output. This output 
can be stored on a hard disk or a floppy disk. 



MEASUREMENT AND COMPUTATION COMPARISON 



10 CENTIMETER THICK WIRE 




Figure 6. Thick wire measurement 
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The impulse response shown in Figure 6 belonged to the signature of the thick copper 
wire having a 90 degree incident wave and a 90 degree reflected wave. The number of 
averages used in the measurement was 2048. The smoothing was C=0.1. There was 
no time drift observed. One data block was used. Other parameters can be found in the 
following table. 



Table 8. 90 DEGREE THICK WIRE PARAMETERS 



Wire length in meters 


0.1 


Wire radius in meters 


0.0026 


Incident angle in degrees 


90 


DG pulse minimum 10% spread in nanoseconds 


0.28 


DG pulse maximum 10% spread in nanoseconds 


0.30 


Time window in nanoseconds 


10.00 


Impact delay in nanoseconds 


2 


Scattering angle in degrees 


90 



Note that the minimum and the maximum 10% spreads determine the respective a, 
and a 2 in the DG waveform of equation 3.9. 

First, the frequency spectrum of the smoothed impulse response was calculated from 
the output and then transformed to the time domain by an FFT [Ref. 18]. The reason 
for calculating the frequency response first is that it is more convenient to work in the 
frequency domain when using Raid's method to extract the smoothed impulse response 
using equation 3.6 

The dotted line in Figure 6 on page 25 is the actual measurement. Because the 
measurement was contaminated by noise, the result did not agree exactly. The contin- 
uous line represents the computed thick wire smoothed impulse response due to the 
same excitation. This computation was performed directly in the time domain using an 
E-field integral equation. Measurements for other types of structures will be illustrated 
in the next chapter. 
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V. EXPERIMENTAL RESULTS 



A. CANONICAL MODEL SCATTERING 

In addition to the 10 cm long thick wire considered in Chapter 4, four additional 
canonical targets were selected: 

• 10 centimeter thin wire (broadside and canted incidence angles) 

• 8.1 centimeter diameter sphere 

• 12.2 centimeter diameter sphere 

• 15.5 centimeter diameter sphere 

The Mie series program, MIE, was used to compute the frequency domain scatter- 
ing transfer function for the various metallic spheres. These transfer functions were then 
used as inputs to the SPRCT program to compute the transient scattering responses due 
to the specified Double-Gaussian (DG) pulse incident plane wave, as specified by 
equation 3.9. The thin wire computations were performed using the TDIG program. 
This program utilizes a numerically evolved time-domain electric field integral equation 
for the induced current distribution on the wire. The scattered field is then calculated 
by appropriate numerical space-time integrations using this induced current. This and 
the sphere scattering programs were written by xM organ. The measured and computed 
DG pulse scattering responses are overlayed in each plot to validate both the fidelity of 
the experimental procedure and the accuracy of the deconvolution algorithm. The 
characteristics of each canonical target scattering case are summarized in Tables 9-13 
on pages 28-33. A noteworthy item is that there is a slight bistatic angle (3 degrees) 
between incident and scattered aspects. This has been incorporated into the computa- 
tions. As can be readily observed, the agreement between measurements and computa- 
tions in Figures 7-11 on pages 28-32 is quite good in all cases. The additive "rippling" 
of the experimental results is due in part to the residual noise and clutter contributions. 
This latter pollutant includes uncancelled multipath interactions of the target and 
canonical sphere scatterers due to imperfect absorption by the anechoic chamber sur- 
faces. Additional differences in the comparisons are due to the approximations which 
are innate to the computations. 
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MEASUREMENT AND COMPUTATION COMPARISON 
12.2 CENTIMETER DIAMETER SPHERE 




TIME SERIES IN NANO SECONDS 



Figure 7. The 12.2 cm diameter sphere validation. 



Table 9. 12.2 CM SPHERE PARAMETERS 



Sphere radius in meters 


0.061 


Distance from sphere in meters 


2.5 


Bistatic angle in degrees 


3 


Scattering plane 


E plane 


Time window in nanoseconds 


20.0 


Narrow 10% pulse width in nanoseconds 


0.28 


Wide 10% pulse width in nanoseconds 


0.30 
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MEASUREMENT AND COMPUTATION COMPARISON 



10 CENTIMETER THIN WIRE 




Figure 8. The 10 cm thin wire validation at 45 degrees 



Table 10. 45 DEGREE THIN WIRE PARAMETERS 



Wire length in meters 


0.1 


Wire radius in meters 


0.00095 


Incident angle in degrees 


45 


DG pulse minimum 10% spread in nanoseconds 


0.28 


DG pulse maximum 10% spread in nanoseconds 


0.30 


Scattering angle in degrees 


135 


No. wire segments 


20 
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MEASUREMENT AND COMPUTATION COMPARISON 



10 CENTIMETER THIN WIRE 




TIME SERIES IN NANO SECONDS 



Figure 9. The 10 cm thin wire validation at 90 degrees 



Table 11. 90 DEGREE THIN WIRE PARAMETERS 



Wire length in meters 


0.1 


Wire radius in meters 


0.00095 


Incident angle in degrees 


90 


DG pulse minimum 10% spread in nanoseconds 


0.28 


DG pulse maximum 10% spread in nanoseconds 


0.30 


No. wire segments 


20 


Scattering angle in degrees 


90 
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MEASUREMENT AND COMPUTATION COMPARISON 
8.1 CENTIMETER DIAMETER SPHERE 




TIME SERIES IN NANO SECONDS 



Figure 10. The 8.1 cm diameter sphere validation 



Table 12. 8.1 CM SPHERE PARAMETERS 



Sphere radius in meters 


0.0405 


Distance from sphere in meters 


2.5 


Bistatic angle in degrees 


3.0 


Time window in nanoseconds 


20.0 


Nyquist frequency in Ghz 


12.775 


Narrow 10% pulse width in nanoseconds 


0.28 


Wide 10% pulse width in nanoseconds 


0.30 
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MEASUREMENT AND COMPUTATION COMPARISON 



15.5 CENTIMETER DIAMETER SPHERE 




TIME SERIES IN NANO SECONDS 



Figure 11. The 15.5 cm diameter sphere validation 



B. SCALE MODEL SCATTERING 

The scattering responses of the wire and sphere targets due to the incident DG in- 
cident pulse can be numerically evaluated using straightforward algorithms. This is not 
the case for more complex scattering objects, such as a highly realistic model of an air- 
craft. Primitive computational algorithms do currently exist for such complex structures 
but the poor accuracy of their results would be insufficient to properly validate the 
measurements being made here. On the other hand, the very good comparisons shown 
in the previous section for both low-damped (wires) and high-damped (spheres) strongly 
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Table 13. 15.5 CM SPHERE PARAMETERS. 



Sphere radius in meters 


0.0775 


Distance from sphere in meters 


2.5 


Bistatic angle in degrees 


3.0 


Scattering plane 


E plane 


Time window in nanoseconds 


20.0 


Nyquist frequency in Ghz 


12.775 


No. frequency points 


257 



support the conjecture that high quality measurements can be obtained for all metallic 
targets having comparable physical dimensions. 

After the validations were made using wires and spheres, extensive measurements 
were made of the transient backscattering from four different silver coated 1/72 scale 
plastic aircraft models. The full-size dimensions of each aircraft, whose identities will 
not be given here, are shown in Tables 14-17 on pages 34-37. 

Measurements were made for various incident aspects in the wing-fuselage plane on 
each model, using horizontal incident polarization: 0 degrees (nose-on); 30 degrees; 90 
degrees (broadside); and 180 degrees (tail-on). In addition, a broadside look-down as- 
pect was considered, with backscattering measured for incident linear polarizations both 
parallel to ("vertical wing") and perpendicular to ("horizontal wing") the fuselage. 
Samples of these numerous measurements are displayed in Figures 12-15 on pages 34-37. 

The complete library of waveforms is being used to support research in demon- 
strating the accuracy and aspect independence of natural resonance radar target iden- 
tification. 

Four scale model aircraft were selected. They are as follows: 

• Scale model 1 

• Scale model 2 

• Scale model 3 

• Scale model 4. 

C. SYSTEM NOISE ESTIMATION 

Noise in the data acquired from the acquisition program was investigated in this 
effort in order to develop a mechanism to estimate the signal-to-noise ratio (SNR). The 
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SCALE AIRCRAFT MODEL 1 



HORIZONTAL WING 





Figure 12. Vertical and horizontal wing scale aircraft model 1 



Table 14. DIMENSIONS OF SCALE AIRCRAFT MODEL 1 



Scale 


1/72 


Length 


36 feet 4 inches 


Wingspan 


31 feet 



34 



SCALE AIRCRAFT MODEL 1 

HEADNQ 000 DECREE WADWO 030 DEGREE 







Figure 13. Scale aircraft model 1 at 000, 030, 090, 180 degree aspect angle 



Table 15. DIMENSIONS OF SCALE AIRCRAFT MODEL 2 



Scale 


1/72 


Overall length 


53 feet 


Extended wingspan 


28 feet 7 inches 



j 
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HEADING 000 DEGREE 

SCALE AIRCRAFT MODEL 1 SCALE AIRCRAFT MODEL 2 







Figure 14. Scale aircraft model comparisons at 000 degree aspect angle 



Table 16. DIMENSIONS OF SCALE AIRCRAFT MODEL 3 



Scale 


1/72 


Length 


56 feet 


Wingspan 


40 feet 8 inches 


Height 


15 feet 4 inches 



36 



HEADING 090 DEGREE 



SCALE AIRCRAFT 1(000. 1 SCALE AJROWT MODQ. 2 






Figure 15. Scale aircraft model comparisons 090 degree aspect angle 



Table 17. DIMENSIONS OF SCALE AIRCRAFT MODEL 4 



Scale 


1/72 


Length 


63.8 feet 


Wingspan 


42.8 feet 


Height 


18.6 feet 
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SNR is a major determining factor for the measurement fidelity. Much of this SNR 
assessment involved the computing some statistics of the waveform. All of these statis- 
tical calculations were performed using GRAFSTAT with APL. 

In the acquisition step, the calibration, the background, the target and a second 
background were measured consecutively. The two backgrounds were measured for the 
purpose of noise estimation. In the following analysis, a 2048 ensemble-average mode 
was used to acquire the data. The acquisition of each waveform took about 9 minutes 
to complete. 

In a separate investigation, the stationarity of the system noise was investigated by 
acquiring seven consecutive background measurements. These backgrounds were then 
subtracted from each other to observe the noise occuring during the consecutive meas- 
urement. For example, the second background and third background were subtracted, 
as is shown in Figure 16 on page 39 along with various statistical parameters in Table 
18 on page 39 mentioned, the statistics of each subtracted waveform was calculated by 
the GRAFSTAT software packgage. The variance of the noise can be calculated from 
the standard deviation and expected power can be approximated by this variance when 
the noise has a small mean. This approximation was reasonable because the squares of 
the mean for the actual measurement data are very small, ( see Figure 16 on page 39 ). 
The computed noise pow'er is plotted versus time in Figure 17 on page 40. As can be 
seen, the measured noise power varied only slightly with time. A linear least-squares lit 
yields the equation 

Y = -93.042 + 0.035727* , (5.1) 

where Y is the expected noise power and X is the time in minutes. This model was 
based on only seven background measurements with only 2048 averages per measure- 
ment. For practical purposes the system noise can be assumed to be fixed in time, de- 
pending only upon the number of ensemble averages used to complete the waveform 
measurement. 

To calculate the SNR, the thick wire measurement was used. The signal plus noise 
waveform is displayed in Figure 18 on page 41 while the computed statistical parameters 
are given in Table 20 on page 41. 

This waveform was time windowed to eliminate the noise outside of the window 
before the statistical parameters were calculated. The signal plus noise power was cal- 
culated by way of the variance using the following equation. 
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EXAMPLE NOISE WAVEFORM 



2ND BACKGROUND MINUS 3RD BACKGROUND 






EATTER PLOT TABLE 


X 


tXXX 


r 


IBCT2-0OT3 


SELECTION 


:ALL 


X LABEL 


:TISE IN NANO SECONDS 


Y LABEL 


: AMPLITUDE IN NANO SEC0MJ6 


ND, OF ELOOTTS 


:912 


CORRELATION XY 


:-0. 34111 


RX CORRELATION 


1-0.14703 T--3.350S 


X MEAN 


:10 


STD. OEVIATION 


:9.7904 


5-PERCENTILE 


to.oazso 


29-PERCEXTILE 


: 4.0902 


MEDIAN 


:10 


75-PERCENTILE 


H9.01 


S3-PERCEKTILE 


:1B.016 


X MIN. 


: 0 0.030139 0.070279 


X MAX. 


:2D 19.861 18.022 


Y MEAN 


J-3.0211E-10 


STD. OEVIATION 


:0. 000023349 


B-PERC ENTILE 


1-0.000030290 


23-PERCOfTILE 


:-7.02E"6 


MEDIAN 


:-7.019E-9 


70-PERCENTILE 


:0. 00002423 


SB-PERCENT I LE 


'.0.000024231 


Y MIN 


j-O. 000069519 -0.000099919 


Y MAX 


: 0.000066731 0.000099461 0, 



Figure 16. Noise waveform statistics 



Table 18. STATISTICS OF NOISE WAVEFORM 



Mean 


-3.021 IE-10 


Std. deviation 


0.000023345 


Median 


-7.01 9E-8 


Maximum 


0.000066731 
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NOISE ANALYSIS 



NOISE POWER OF SIX TRIAL 




SCATTER PLOT TABLE 



X 


TIME 


Y 


10x10»PCWERS 


SELECTION 


ALL 


X LABEL 


TIME IN MINUTE 


Y LABEL 


DB 


NO. OF ELEMENTS 


6 


CORRELATION XY 


0.60163 


RK CORRELATION 


0.63775 T-1.658 


X MEAN 


31.5 


STD. DEVIATION 


16.637 


5-PERCENTILE 


9 


25-PERCENTILE 


16 


MEDIAN 


31.5 


75-PERCENTILE 


45 


95-PERCENTILE 


54 


X MIN. 


9 18 27 


X MAX. 


54 45 36 


Y MEAN 


“91.917 


STD. DEVIATION 


0.999B6 


5-PERCENTILE 


“92.76 


25-PERCENTILE 


“92.716 


MEDIAN 


“92.343 


75-PERCENTILE 


“90.669 


95-PERCENTILE 


“90.669 


Y MIN 


“92.78 “92.718 “92.638 


Y MAX 


“90.669 “90.669 “92.049 



Figure 17. Noise power 



Table 19. STATISTICS OF NOISE MODEL 



Mean 


-91.917 


Std. deviation 


0.99986 


Median 


-92.343 


Maximum 


-90.669 



40 



THICK WIRE SIGNAL PLUS NOISE 




SCATTER PLOT TABLE 
X :XXX 

Y :TOTALx(rmiCX-TTHIBCT) 

SELECTION :ALL 

X LABEL :TIME IN NANO SECONDS 

Y LABEL : AMPLITUDE IN NANO SEOWCS 

NO. OF ELEMENTS :512 

CORRELATION XY : 0.0009901 
RK CORRELATION :0. 00*1907 T^>. 13901 
X LEAN *10 

STD. DEVIATION : 5. 7904 
5-PERCENTILE : 0.90239 
25-PERCENTILE : 4. 9902 
MEDIAN : 10 

75-PERCENTILE : 15.01 
95-PERCENTILE : 19. 010 
X MIN. :0 0.039139 0.070270 

X MAX. :20 19.901 19.922 

Y LON :-0.2B91E-0 

STD. DEVIATION : 0.00033 62 
5-PERCENTILE : **0.00058722 
25-PERCENTILE :0 

MEDIAN :0 

75-PERCENTILE :0 
95-PERCENTILE : 0.00069710 

Y MIN j-O. 0014622 “0.001300 “0.0012747 

Y MAX :0. 0014120 0.001319 0.0012253 



Figure 18. Signal plus noise waveform statistics 



Table 20. STATISTICS OF THICK WIRE SIGNAL PLUS NOISE 



WAVEFORM 



Mean 


-6.299 IE-6 


Std. deviation 


0.0003362 


Median 


0 


Maximum 


0.0014128 
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E(X 2 ) = Variance + E 2 (X) 



(5.2) 



where E(X) is mean of the noise and E{X 2 ) is the mean square value of the noise (noise 
power). Note that the second term on the right side of equation 5.2 was close to zero, 
so it can be neglected. As a result, the expected power is equal to the variance. The 
signal plus noise power was equal to 0.0003382 2 or -69.41 dB. The noise expected power 
was calculated by equation 5.1 and is given by -92.64 dB for the average time of 10 
minutes. Since the noise power was so small compared to the signal plus noise power, 
the subtraction between the signal plus noise and the noise was not necessary. The 
noise power when time windowed was reduced by a factor of 200/512, or 4.08 dB. The 
result for the SNR is 27.31 dB in the case of the thick wire measurement. Note that this 
SNR is about 10 dB better than that provided for the same target using the Tektronix 
equipment, as carried out by McDaniel [Ref. 11]. 
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VI. CONCLUSIONS 



A. SUMMARY 

This thesis describes the theory and implementation of a newly updated transient 
scattering measurement facility. An overview of the historical precedents and moti- 
vations for this development was given in Chapter I. This was followed in Chapter II 
by a detailed description of the new equipment being used as well as the physical struc- 
ture of the anechoic chamber. 

Building upon this overview of the hardware aspects, the conceptual and math- 
ematical modeling of the scattering range is undertaken in Chapter III. Using the 
physical model of the radiating and receiving antennas in the anechoic chamber, the 
various electromagnetic interactions between the scattering target and the chamber are 
represented by a comprehensive linear system model. The topology of the various fre- 
quency domain transfer functions in this model is then used to justify the measurement 
and pre-processing steps required to estimate the smoothed scattering impulse response, 
due to a specified double-Gaussian incident plane wave. There are three measurements 
required for each target impulse response estimation. Aside from the actual target 
measurement, a background measurement is needed to subtract out the directly received 
signal from the transmitting antenna, in addition to that scattered from the various 
chamber surfaces. A third measurement, for a metallic sphere, is needed to allow for 
compensation of the transfer function characteristics of the transmitting and receiving 
antennas in the chamber. By effecting a deconvolution comparison with a highly accu- 
rate numerical calculation of the transient scattering response of the sphere target, this 
final measurement also permits the estimation of the target's response to the same 
excitation as is used in the sphere computation. 

Chapter IV considers the details of the actual software that performs the signal ac- 
quisition, signal averaging and deconvolution post-processing of the measured 
waveforms. The acquisition and signal averaging software is a blend of a commercial 
IEEE-488 bus controller (GURU-II) and a custom driver, designed as part of this thesis 
effort. The deconvolution algorithm provides, first of all, graphically displayed iterative 
background subtractions from both the target and sphere measurements to allow com- 
pensation for time drifts. Secondly, the final deconvolution is performed via a variant 
on Riad's method, as discussed in Chapter III. 
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The validations shown in Chapter V speak for themselves as to the expected meas- 
urement quality obtained using this new transient scattering range. Both low-Q and 
high-Q type targets are considered by the respective comparisons of sphere and thin-wire 
scattering measurements with computed smoothed impulse responses. Further illus- 
trations of the type of measurements to be performed with this range are shown in the 
scale model aircraft backscattering signatures. A simple software modification will allow 
the display of very broadband frequency domain radar cross sections (RCS) for scale 
model targets. 

To quantitatively investigate the fidelity of the measurement system, a noise esti- 
mation of the new implementation was investigated and was modelled. The particular 
signal-to-noise ratio (SNR) of a thick wire scattering measurement was calculated. A 
27 dB SNR of this case demonstrated the relatively high fidelity of the measurement 
system for the case small RCS targets, including silver coated 1/72 scale model aircraft. 

B. FUTURE CONSIDERATIONS 

Even though the demonstrated system performance was excellent, there is no such 
thing as too high an SNR in scattering measurements. In this regard, the power- 
bandwidth product has the potential to be enhanced by future efforts. 

The 1 to 7 GHz passband of the scattering system was limited by the transmitting 
GaAS power amplifier. On the other hand, the bandwidth of the HP DPO system can 
go as high as 20 Ghz, while the current antennas have a 1 to 12 GHz passband. Future 
use of broader bandwidth antennas and power amplifier, having a rated RMS output 
power of 2 watts or more, would increase the overall system performance. More accu- 
rate smoothed impulse responses could be achieved, as well as wider bandwidth RCS 
measurements. 

Because of the limited scope of this thesis, the noise power estimation was only an 
approximation. A more comprehensive noise analysis should be carried out in order to 
quantify and minimize the noise power sources. 

The final consideration is the future investigation of alternate mechanisms for pro- 
viding broadband scattering measurements. One such mature technique is through a 
stepped frequency continuous wave system, employing a network analyzer front-end as 
the coherent receiver, perhaps preceded by a low-noise amplifier to reduce the overall 
system noise figure. A somewhat radical technique, one that is currently being investi- 
gated at NPS, is to use a very broadband amplified noise source as the transmitting 
source, while employing the HP DPO as a two-channel correlation receiver. This 
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method is akin to a dynamic matched filter, offering enhanced bandwidth, at relatively 
low cost vis-a-vis the stepped frequency setup. Additionally, the potential exists for both 
tactical and strategic employment of such random source systems, either as active 
"pseudo-jammers" or as totally passive radars relying upon either natural or man-made 
random signal sources for the illumination of targets. 
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APPENDIX A. THE IMPULSE RESPONSE DECONVOLUTION 



The following computer program calculates the Impulse Response from the meas 
urement data as described in chapter 4. The program is written in RM/FORTRAN. 

The subroutine PL08 and Pll require the GRAFMATIC routine library. The pro 
cedure to compile, link and run this source code refers to the RM/FORTRAN User' 
Guide Version 2.11 (DOS). 

£************************************************* 

C THIS IS THE IMPULSE RESPONSE DECONVOLUTION * 

C WRITTEN BY LT. JG SOONPUEN SOMAPEE Feb 88 * 

C ADAPTATION OF ORIGINAL TEKTRONIX PROGRAM * 

C BY PROF. M. A. MORGAN Circa 1983-85 * 

Q*** * * *** **** ** * * * * * * ** ***** ** ** *** ** ** ** * * ** ** *** 
Q************************************************* 



c 


WHERE: 




* 


c 


DATE 


IS THE DATE OF MEASUREMENT 


* 


c 


DATAB 


IS THE NUMBER OF DATA BLOCKS 


* 


c 


FI 


IS THE DUMMY STRING 


* 


c 


ID 


IS THE IDENTIFICATION STRING 


* 


c 


K 


IS THE NUMBER OF POINTS 


* 


c 


L9 


IS THE RELATIVE SMOOTHING 


* 


c 


MAG 


IS THE MAGNITUDE OF Xo*Ht(f) 


* 


c 


NAME4 


IS THE NAME OF THE TARGET DATA 


* 


c 


N 


IS THE NUMBER OF POINTS 


* 


c 


NUMPT 


IS THE SAME AS K 


* 


c 


NEXPT 


IS INDICATEES POWER-OF-2 EXPONENT 


* 


c 




WILL BE GENERATE BY PREFFT 


* 


c 




SET TO -1 INDICATES ERROR 


* 


c 


R 


IS THE REAL PART OF THE INVERSE 


* 


c 




FOURIER TRANSFORM OF Xo*Ht(f) 


* 


c 


SUBAV 


IS THE NUMBER OF SUB -AVERAGES 


* 


c 


T 


IS THE TIME INTERVAL 


* 


c 


TIME 


IS THE TIME INTEVAL 


* 


c 


TITLE 


IS THE TITLE NAME FOR PLOTTING 


* 


c 


TUT 


IS THE FREQUENCY DOMAIN OF Xo*H(f) 


)* 


c 


TYPE 


IS THE WAVEFORM TYPE 


* 


c 


W 


IS THE COMPLEX EXPONENTIAL ARRAY 


* 


c 




GENERATED BY PREFFT SUBROUTINE 


* 


c 


Xo*Ht(f) 


IS THE GAUSSIAN WEIGHTED COSINE 


* 


c 


xs 


IS THE STRING FOR PLOTTING ROUTINE* 


c 




FOR THE HEADER 


* 


c 




DRIVER 


* 


c 


Y4 


IS THE CALIBRATION SPHERE 


* 


c 

n 




SUBTRACTED WAVEFORM (FREQUENCY) 


* 

* 


\J 

c 


Y5 


IS THE TARGET SUBTRACTED WAVEFORM 


* 


c 

r 




FREQUENCY DOMAIN 


* 

* 


\J 

c 


YSC 


IS THE COMPUTED SCATTERED ELECTRIC* 
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c 

c 


Y1 


FIELD FOR THE CANONICAL TARGET 
IS THE CALIBRATION MEASUREMENT 


* 

* 


c 


Y2 


IS 


THE 


BACKGROUND MEASUREMENT 


* 


c 


Y3 


IS 


THE 


TARGET MEASUREMENT 


* 


c 


YS 


IS 


THE 


SAME AS XS 


* 


c 


YN 


IS 


THE 


YES OR NO STRING 


* 


c 


XMI 1 


IS 


THE 


MINIMUM VALUE 


* 


c 


XMA 1 


IS 


THE 


MAXIMUM VALUE 


* 



CHARACTER*64 TITLE , XS , ID , YS , DATE ,TYPE , NAME 
CHARACTERS TYPE1 ,TYPE2,TYPE3,F1 
CHARACTER* 16 NAME 1 , NAME2 , NAME3 
CHARACTERS ID1,ID2,ID3 
CHARACTER* 1 YNS,YNS1 

COMPLEX YSC( 1025) , Y4( 1025) ,Y5( 1025) ,W( 1025) 
REAL T,T2 ,TIME 

REAL Yl( 1024), Y2( 1024), Y3( 1024) 

INTEGER N , SUB AV , DATAB , NUMPT 

Q************************************************* 

YNS='N’ 

Q * * * * * * * * * * ** * * * * -V * ** * * * * * * ** ****** ** ************* 



C READ DATA FOR COMPUTED SPHERE THEN 
C CALLING SUBROUTINE FFF TO DRAW THE TIME 
C SERIES AND .SPECTRAL PLOT 

Q * * V: * * * * “V * * * * * -V * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 

WR I TE ( * 5 * ) f **************************************************** 



*) 


T * 










* 


*) 


1 * 


DECONVOLUTION PROGRAM 


BY DR. M. A. MORGAN 


* 


*) 

*) 

*) 

*) 

*) 


f * 


ADAPTED 


TO IBM BY LT. JG. 


SOONPUEN SOMAPEE 


* 


1 * 






FORMAT OF DATA FILE: 


* 


f * 


1. 


READ 


ID 


(FORMAT CHARACTER* 64) 


* 


*) 


! * 


2 . 


READ 


TYPE 


(FORMAT CHARACTER* 64) 


* 


*) 


t y- 


3. 


READ 


DATE 


(FORMAT CHARACTER* 64) 


* 


*) 


T * 


4. 


READ 


FI 


(FORMAT CHARACTER*64) 


* 


*) 


1 * 


5. 


READ 


SUB -AVERAGE 


(FORMAT INTEGER) 


* 


*) 


T * 


6 . 


READ 


TIME WINDOW 


(FORMAT REAL(F5. 2)) 


* 


*) 


T * 


7. 


READ 


DATA BLOCK 


(FORMAT INTEGER) 


* 


*) 


T * 


8 . 


READ 


NUMBER OF POINTS 


(FORMAT INTEGER) 


* 



WRITE(*, 
WRITE(* , 
WRITE (*. 



WR I TE ( * * ) * **************-V**-V*******-!Wr************Vr**V 
WRITE ( ,v * ,v ) ’FILE NAME FOR COMPUTED CALIBRATION SPHERE: 
READ(*,173) NAME 

Q * * * * * * * * * * * * * * * * * * * * * ********** * * ***** * * * ********** * * * * * * * * 



?******** 



c 

C YSC = COMPUTED SPHERE COMPLEX FREQUENCY DOMAIN DATA 
C 

C NUMPT = NUMBER OF POINTS IN THE OUTPUT OF THIS ROUTINE 
C 

C TIME = OUTPUT TIME INTERVAL 
C 

c*********************************************************** 
CALL FFF 1( NAME, YSC, W,NEXPT, TIME) 
c*********************************************************** 
9110 CONTINUE 

c*********************************************************** 
C READ THE DATA FILES FOR CALIBRATE AND BACKGROUND 
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C THEN OVERLAY THESE ONTO THE SAME AXIS. 

C PLOT77 ALLOWS THE OPERATOR TO SHIFT ONE WAVEFORM 
C PRIOR TO SUBTRACTING THE BACKGROUND. 

C 

C THE RESULT FROM SUBTRACTION IS TRANSFORMED TO 
C FREQUENCY DOMAIN AND SAVED TO THE COMPLEX ARRAY Y4 
C 

C NAME1 = NAME OF THE FILE TO SENT TO PLOT77 ROUTINE 
C (CALIBRATE DATA) 

C NAME 2 = NAME OF THE BACKGROUND DATA 

C ID1 = IDENTIFICATION CHR$ FOR CALIBRATE 

C Y4 COMPLEX ARRAY OUTPUT (CALIBRATE ) 

Qkkkkkk k k * kkkk k kk kkkkkk kk kk kkk kk kkkk kkkkkk kkkkkkkkkkkkkkk kkk 

WR I TE ( * ) 1 * ’V k kkk kkkk k Vc-iY kk * Vc * kkkkkkk kkk kkkkkkkk kkkk * 1 

WRITE(*’*) ’FILE NAME FOR CALIBRATION SPHERE: ’ 

READ(*,173) NAME1 

CALL I NPP ( NAME 1 , 1 D 1 , TYPE 1 , DATE ,F1,SUBAV,T2, D ATAB , NUMPT , Y 1 ) 
WR I TE ( * * ) 1 * Vr * * *>WnY k 'kick ft kkk Vr* *>V Vr ***** VbY kkkkkk Ve *?Y kkkk *>Wc ! 

WRITE(*|*) 'FILE NAME FOR SPHERE BACKGROUND WAVEFORM: ’ 
READ( ,V , 173) NAME 2 

C * * *** * irk k k k irk kkk k k kk kkk kk kkkkkk k kkk k kk kk kk kk kkkk 

CALL INPP( NAME 2 , ID2 , TYPE 2 , DATE , F 1 , SUB AV , T2 , DATAB , NUMPT , Y2 ) 
CALL PL077(Y1, Y2,T2, NUMPT, ID1) 

Qkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 



9112 CONTINUE 

WR I TE ( * y k ) 1 k k k k k k k k k k k k k k k k kk kkkk k kkkk kkkkkkk kkkk kkk 
WRITE (*V v ) 'FILE NAME FOR TARGET WAVEFORM: ’ 

READ(* ,173) NAME 3 



Q k kk k k k k k k k k k k k kk k k kk kk kk k k k k k k k k k kk kkk k k kk k k k kkk 



i 



CALL I NPP( NAME 3 , ID3 , TYPE 3 , DATE , F 1 , SUB AV , T2 , DATAB , NUMPT , Y3 ) 

0 k kk kkk kkkk k kkkkk kkkkkkkkk k kkkk kkk kkk k kk kk kk kk k kk 



WRITE(*,*) 

WRITE (*,*) 'USE PREVIOUS BACKGROUND FOR THE TARGET ? (Y/N): 
READC*, 173) YNS1 

IF(YNS1.EQ. 'Y’.OR. YNS1.EQ. V) GO TO 9222 

WRITE ( * * ) 1 kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 1 

WRITERS*) 'FILE NAME FOR TARGET BACKGROUND WAVEFORM: ’ 

READ(*, 173) NAME 2 

0 k k kk k kk k k k k k k k k k k k k k kk k k k k k k kk kk k kk kk k kk k kkkk k k k 

CALL INPP( NAME 2 , ID2 , TYPE 2 , DATE , FI , SUBAV , T2 , DATAB , NUMPT , Y2 ) 

9222 CONTINUE 
C 

C REPEATING SHIFT AND SUBTRACT OPERATIONS ON TARGET AND BACKGROUND 
C NAME1= TARGET DATA FILE 

C NAME 2= BACKGROUND DATA FILE 

C OUTPUT 

C ID1 = IDENTIFICATION OF THE TARGET 

C Y5 = COMPLEX ARRAY OF FREQUENCY DOMAIN TARGET DATA 

Qkkkkkk k k kkkkkkkk k kkkkkk k k k k k k k k k k k k kkkk k kkkkkkkkk 



CALL PL077( Y3 , Y2 ,T2 , NUMPT, ID3) 

DO 987 1=1, NUMPT 

Y4( I)=CMPLX(Y1( I) ,0. 0) 

Y5( I )=CMPLX(Y3( I) , 0. 0) 

987 CONTINUE 

MODE=0 

CALL FFT( NUMPT , MODE , TIME , NEXPT , W , Y4 ) 
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CALL FFT( NUMPT , MODE , TIME , NEXPT , W , Y5 ) 

WRITE ( * * ) ' * * * * * * * * * * 'V * * * * * * * * * ** * * * iV * * * * * * iV * * * * iV * * * 1 
CALL DECON3 ( NUMPT , YSC , Y4 , Y5 , ID3 , TYPE 3 , DATE , F 1 , SUB AV , T2 , DATAB ) 

WRITE(*’*) 'WOULD YOU LIKE TO PROCESS ANOTHER TARGET ? (Y/N): ' 
READ(*,173) YNS 

IF(YNS. EQ. ' Y' .OR. YNS.EQ. 'y' ) GOTO 9112 
VR I TE ( * 5 * ) f 1 

WRITE (*,*) ' * eeeee eeeeee e *' 

WRITE(*,*) '* eee eee eeee *' 

WRITE(*,*) ' eeeeeeeeeeeeeeeeeeeeeeeeeeeee ,v ' 

WRITE(*,*) ' eeeeeeeeeeeeeeeeeeeeeeee *' 

WRITE(*,*) 1 eeeeeeeeeeeeeeeeeeeeee *' 

WRITE( * * ) 1 ******* * * * * * * * * * * * * * * * * * * * * * * * * * * 1 
WRITE ( * , * ) * ********************************* ' 

173 FORMAT(A) 

STOP 

END 



SUBROUTINE INPP( NAME, ID, TYPE, DATE, FI , SUB AV,T2, DATAB, NUMPT, Y) 
q*********************** ************************************** ******** 



C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

0 ****^ 



INPUT 
OUTPUT 
NAME = 
ID 

TYPE = 
DATE = 
INTEGER 
SUBAV = 
DATAB = 
FI 

NUMPT = 
REAL 
T2 
Y 



NAME 

ID , TYPE , DATE , F 1 , SUBAV , T2 , DATAB , NUMPT , Y 
CHR$ NAME OF THE DATA FILE 
CHR$ TARGET IDENTIFICATION 
CHR$ WAVE FORM IDENTIFICATION 
: CHR$ DATE 

NUMBER OF SUB -AVERAGES USED TO ACQUIRE THE DATA 

DATA BLOCK 

DUMMY 

NUMBER OF POINTS 

TIME WINDOW 
Y ARRAY OUTPUT 



A 4 \ 4 Z 4\1\ 4\ 4\ 4\7Z 4\ /V 4\ 4\ 4 Z 






;V it it iz Vr it * * it V r iz * V? * * it Vr V? ** iz Vc * it Vr * Vc Vr iz iz Vc Vr V: it V? * Vc * * * iz V? i: iz 



CHARACTER* 16 NAME 
CHARACTER* 6 4 ID, FI 
CHARACTER*64 TYPE, DATE 
REAL* 4 Y( 1500) 

INTEGER SUBAV, DATAB, NUMPT, K 
PI=3. 141592654 



OPEN (UNIT=10, FILE =NAME,STATUS=' OLD’ ) 
C READ THE HEADER 

READ( 10,5) ID 
READ( 10,5) TYPE 
READ( 10,5) DATE 
READ( 10,5) FI 
READ( 10 , 2030) SUBAV 
READ( 10,2010) T2 
READ( 10,2030) DATAB 
READ( 10,2030) NUMPT 
10 READ( 10 , 2000 ,END=20) Y(K) 
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K=K+1 


17 


GOTO 10 


20 


CLOSE(IO) 


2000 


FORMATCE12. 6) 


2010 


FORMAT(F5. 2) 


2030 


FORMAT( 15) 


5 


FORMAT(A) 



K=K- 1 

IF(K. EQ. NUMPT) GOTO 3000 

WRITE(*,*) ’CHECK THE DATA - # POINTS MISMATCH' 
STOP 

3000 RETURN 
END 

SUBROUTINE FFF1( NAME, YY,W,NEXPT, TIME) 



C 


WHERE 




* 


C 


DATE 


IS THE DATE OF MEASUREMENT 


* 


c 


DATAB 


IS THE NUMBER OF DATA BLOCKS 


* 


c 


FI 


IS THE DUMMY STRING 


* 


c 


ID 


IS THE IDENTIFICATION STRING 


* 


c 


K 


IS THE NUMBER OF POINT 


* 


c 


L9 


IS THE RELATIVE SMOOTHING 


* 


c 


MODE 


IS EQUAL TO 0 MEAN FORWARD 


* 


c 




FOURIER TRANSFORM 




c 


MAG 


IS THE MAGNITUDE OF Xo*Ht(f) 


* 


c 


NAME 4 


IS THE NAME OF THE TARGET DATA 


* 


c 


N 


IS THE NUMBER OF POINTS 


* 


c 


NUMPT 


IS THE SAME AS K 


* 


c 


NEXPT 


IS THE POWER -OF -2 EXPONENT 


* 


c 




WILL BE GENERATED BY PREFFT 


* 


c 




SET TO -1 INDICATES ERROR 


* 


c 


R 


IS THE REAL PART OF THE INVERSE 


* 


c 




FOURIER TRANSFORM OF Xo*Ht(f) 


* 


c 


SUBAV 


IS THE NUMBER OF SUB AVERAGE 


* 


c 


T 


IS THE TIME INTERVAL 


* 


c 


TIME 


IS THE TIME INTEVAL 


* 


c 


TITLE 


IS THE TITLE NAME FOR PLOTTING 




c 


TUT 


IS THE FREQUENCY DOMAIN OF Xo*H(f)* 


c 


TYPE 


IS THE WAVEFORM TYPE 


* 


c 


W 


IS THE COMPLEX EXPONENTIAL ARRAY 


* 


c 




GENERATE BY PREFFT SUBROUTINE 


* 


c 


Xo*Ht( f ) 


IS THE GAUSSIAN WEIGHTED COSINE 


* 


c 


xs 


IS THE STRING FOR PLOTTING ROUTINE* 


c 




FOR THE HEADER 


* 


c 




DRIVER 


* 


c 


Y4 


IS THE CALIBRATION SPHERE 


* 


c 

p 




SUBTRACTED WAVEFORM (FREQUENCY) 


* 

* 


c 


Y5 


IS THE TARGET SUBTRACTED WAVEFORM 




c 

p 




FREQUENCY DOMAIN 


* 

* 


c 


YSC 


IS THE COMPUTED SCATTERED ELECTRIC* 


c 




FIELD FOR THE CANONICAL TARGET 


* 


c 


Y1 


IS THE CALIBRATION MEASUREMENT 


* 


c 


Y2 


IS THE BACKGROUND MEASUREMENT 
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c 


Y3 


IS THE TARGET MEASURED 


* 


c 


YS 


IS THE SAME AS XS 


* 


c 


YN 


IS THE YES OR NO STRING 


* 


c 




METER 


* 


c 


XMI 1 


IS THE MINIMUM VALUE 


* 


c 


XMA1 


IS THE MAXIMUM VALUE 


* 



CHARACTER* 16 NAME 

CHARACTER* 6 4 I D , TYPE , DATE , TITLE , XS , YS , F 1 
COMPLEX W( 1025), YY( 1025) 

REAL T 

REAL Y( 1024) , TIME ,T2 ,XMINS ,MAG( 1024) ,XMAG 
INTEGER NUMPT,DATAB, SUB AV, MODE, N 

Qick'k'k'k •k Vc tVVcjV * * * * * * * * * ?V * * * ** V«V ** *>V ** * * -jV?V * Vc * VoV * tV-jV^V * 

CALL INPP( NAME, ID, TYPE, DATE, FI, SUB AV,T2,DATAB,NUMPT,Y) 

TIME=T2/(FL0AT(NUMPT-l))*lE-9 

XMINS=0. 0 

N=NUMPT 

CALL PRETIM( TITLE , XS , YS ) 

CALL P 11( TITLE ,N,XMINS ,T2 ,Y,XS, ID, YS) 

C SET MODE EQUAJ TO 0 TO FIND FOURIER SERIES 
MODE=0 

CALL PREFFT ( N , MODE , NEXPT , W ) 

DO 1234 I =1 , NUMPT 

YY( I )=CMPLX( Y( I) ,0. 0 ) 

1234 CONTINUE 

T=TIME 

CALL FFT (N, MODE, T, NEXPT, W,YY) 

C CHANGING THE OUTPUT YY TO MAGNITUDE AND PHASE 
N=NUMPT 

CALL PREFRE( TITLE, XS,YS) 

CALL MAGN I TUD ( YY , N , MAG ) 

XMAG=1. O/TIME/2. 0*lE-9 

CALL P 11( TITLE, N, XMINS ,XMAG, MAG, XS , ID , YS) 

RETURN 

END 



C* 

c 

c 

c 

c 

c 

c 






* 



FFT PROGRAM FROM THE DIGITAL SPECTRAL ANALYSIS * 
WITH APLICATIONS PRENTICE -HALL, INC. . * 

ENGLEWOOD CLIFF, NEW JERSEY 07632 1987 * 

S. LAWRENCE MARPLE * 

* 






SUBROUTINE PREFFT (N,MODE ,NEXP,W) 



C 

C 

C 

C 

C 

C 

C 

C 

C 

C 



Input Parameters: 

N - Number of data samples to be processed (integer-must be a 
power of two) 

MODE - Set to 0 for discrete-time Fourier series (Eq. 2. C. 1) or 
1 for inverse (Eq. 2. C. 2) 

Output Parameters: 
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NEXP - Indicates power-of-2 exponent such that N=2 ,V,V NEXP . 

Will be set to -1 to indicate error condition if N 
is not a power of 2 (this integer used by sub. FFT) 

W - Complex exponential array 

Notes: 

External array W must be dimensioned . GE. N by calling program. 

COMPLEX W( 1024) ,C1,C2 
NEXP=1 

5 NT=2**NEXP 

IF (NT .GE. N) GO TO 10 

NEXP=NEXP+1 

GO TO 5 

10 IF (NT .EQ. N) GO TO 15 

NEXP=-1 
RETURN 

15 S=8. *ATAN( 1. ) /FLOAT (NT) 

C1=CMPLX(C0S(S),-SIN(S)) 

IF (MODE . NE. 0) C1=C0NJG(C1) 

C2=( 1. ,0. ) 

DO 20 K=1,NT 
W(K)=C2 

20 C2=C2*C1 

RETURN 
END 



Q Vc •> V -jV •> V Vr Vc Vc Vr * Vr Vr * Vr “>V Vc •> V Vc Vc Vc * -jV Vr Vr Vc *>Y Vc -)Y * Vc •>> *>V Vc * * ->Y - Y * Vc ■> V Vc ->Y Vc Vc *>Y Vc 



SUBROUTINE FFT (N, MODE, T, NEXP, W,X) 

C 

C Input Parameters: 

C 

C N, MODE , NEXP ,W - See parameter list for subroutine PREFFT 

C T - Sample interval in seconds 

C X - Array of N complex data samples, X(l) to X(N) 

C 

C Output Parameters: 

C 

C X - N complex transform values replace original data samples 

C indexed from k=l to k=N, representing the frequencies 

C (k-l)/NT hertz 

C 

C Notes: 

C 

C External array X must be dimensioned . GE. N by calling program. 
C 

COMPLEX X( 1024) ,W( 1024), Cl, C2 

MM=1 

LL=N 

DO 70 K=1 ,NEXP 
NN=LL/2 
JJ=MM+1 

DO 40 1=1 ,N ,LL 
KK=I+NN 
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C1=X( I)+X(KK) 

X(KK)=X( I) -X(KK) 

40 X(I)=C1 

IF (NN .EQ. 1) GO TO 70 
DO 60 J=2 ,NN 
C2=W( JJ) 

DO 50 I=J,N,LL 
KK=I+NN 
C1=X(I)+X(KK) 

X(KK)=(X( I) -X(KK) )*C2 
50 X(I)=C1 

60 JJ=JJ+MM 

LL=NN 
MM=MM*2 

70 CONTINUE 

NV2=N/2 
NM1=N- 1 
J=1 

DO 90 1=1, NM1 

IF (I . GE. J) GO TO 80 

C1=X(J) 

X(J)=X(I) 

X(I)=C1 
80 K=NV2 

85 IF (K .GE. J) GO TO 90 

J=J-K 
K=K/2 
GO TO 85 
90 J=J+K 

IF (MODE .EQ. 0) S=T 
IF (MODE .NE. 0) S=l. /(T*FLOAT(N) ) 
DO 100 1=1, N 
100 X(I)=X(I)*S 

RETURN 
END 



Q Vc 'k Vc *>V * jc Vc Vc *>V Vc Vr * Vr Vc jc *>V *>V * *>V *>V "/r>V Vr Vc tV -k •k *>V Vc * ?Wr Vr * *>V 

c 

SUBROUTINE MAGNITUD ( Y , NUMPT , MAG ) 

C 

Q ■> V *>V *>V * Vr Vc Vc *>V ?V Vc *>V Vc Vc * Vr Vr * •jV * * * Vc Wc ■) V Vo Wr *>V * Vc ? V Vc Vc * * tV tV *>V * V? * *>V *>V * 

COMPLEX Y( 1024) ,H 
REAL*4 MAG( 1024) 

INTEGER N,K,NP, NUMPT 
DO 10 1=1 , NUMPT / 2 
H=Y( I) 

MAG( I )=CABS(H) 

10 CONTINUE 

NUMPT=NUMPT / 2 

RETURN 

END 



QVfycyc*****v< , yc**Vf***7V****Vf , 3V^c , *yf **************’/?'/?******* 



SUBROUTINE PRETIM( TITLE ,XS ,YS) 
CHARACTER*64 TITLE, XS,YS 
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TITLE=' TIME SERIES PLOT' 

XS =' X-AXIS TIME IN NS ' 

YS =' Y-AXIS MAGNITUDE IN VOLT' 

RETURN 

END 



0 ************************************************* 



SUBROUTINE PREFRE ( TITLE ,XS ,YS) 
CHARACTERS TITLE, XS,YS 
TITLE=' SPECTRAL PLOT' 

XS=' X-AXIS FREQUENCY IN GHZ' 
YS=' Y-AXIS MAGNITUDE IN VOLT' 
RETURN 
END 



SUBROUTINE IFFF(N,T,X) 

0 ** **** * *** ********* * ** * * ***** ***************************** 
COMPLEX W( 1025 ),X( 1025) 

REAL T 
INTEGER MODE 
M0DE=1 

CALL PREFFT(N,MODE ,NEXP,W) 

CALL FFT ( N , MODE , T , NEXP , W , X) 

RETURN 

END 



Q********************************************************************* 



SUBROUTINE DEC ON 3 (K, YSC , Y4,Y5 , ID, TYPE , DATE , F 1 , SUB AV , T2 , D ATAB ) 
0 * * ***** ********* ******* ** **** *** ***** *** * *********** 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



COMPUTING THE TARGET RESPONSE TO THE INCIDENT * 
WAVESHAPE FOR THE COMPUTED SPHERE. USING * 
RIAD’S METHOD OF OPTIMAL DECONVOLUTION. * 

* 



Xo*Ht(f)=Y5*C0NJ(Y4)*Ysc/(Y4*C0NJ(Y4)+L) 

WHERE 

Xo( f ) IS THE INCIDENT WAVEFORM 

Ht(f) IS THE TARGET TRANSFER FUNCTION 

Y4 IS THE CALIBRATION SPHERE 

SUBTRACTED WAVEFORM (FREQUENCY) 

Y5 IS THE TARGET SUBTRACTED WAVEFORM 

FREQUENCY DOMAIN 

YSC IS THE COMPUTED SCATTERED ELECTRIC 

FIELD FOR THE CANONICAL TARGET 
Y1 IS THE CALIBRATION MEASUREMENT 

Y2 IS THE BACKGROUND MEASUREMENT 

Y3 IS T?'E TARGET MEASURED- 

K IS THE NUMBER OF POINT 

NUMPT IS THE SAME AS K 

R IS THE REAL PART OF THE INVERSE 



* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

->v 

* 

* 
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c 




FOURIER TRANSFORM OF Xo*Ht(f) 


* 


c 


MAG 


IS 


THE 


MAGNITUDE OF Xo*Ht(f) 


* 


c 


ID 


IS 


THE 


IDENTIFICATION STRING 


* 


c 


SUBAV 


IS 


THE 


NUMBER OF SUB AVERAGE 


* 


c 


DATAB 


IS 


THE 


NUMBER OF DATA BLOCK 


* 


c 


FI 


IS 


THE 


DUMMY STRING 


* 


c 


NAME 4 


IS 


THE 


NAME OF THE TARGET DATA 


* 


c 


XS 


IS 


THE 


STRING FOR PLOTTING ROUTINE 


* 


c 




FOR THE HEADER 


* 


c 


YS 


IS 


THE 


SAME AS XS 


* 


c 


TITLE 


IS 


THE 


TITLE NAME FOR PLOTTING 


* 


c 


YN 


IS 


THE 


YES OR NO STRING 


* 


c 


L9 


IS 


THE 


SMOOTHING PARAMETER 


* 


c 


TUT 


IS 


THE 


FREQUENCY DOMAIN OF Xo*H(f) 


* 


c 


TYPE 


IS 


THE 


WAVEFORM TYPE 


* 


c 


XMI1 


IS 


THE 


MINIMUM VALUE 


* 


c 


XMA1 


IS 


THE 


MAXIMUM VALUE 


* 


c 


DATE 


IS 


THE 


DATE OF MEASUREMENT 


* 



Q ******* * *** *** *** ***** *** * ** * * * * * * * * * **** *** *** ***** 



COMPLEX 

COMPLEX 

REAL 

REAL 

INTEGER 

CHARACTER* 16 

CHARACTER*64 

CHARACTER*! 



Y5( 1025) , Y4( 1025) ,YSC( 1024) ,TUT( 1024) , SUM 
LI ,L2 ,L3 

R( 1024), L9,XMA1,XMI1,MAG( 1024) 

TIME,T,T2 

K , NUMPT , SUB AV , DATAB 
NAME 4 

ID, YS,XS, TITLE, TYPE, DATE, FI 
YN 

c************************************************ 
T=T2/FLOAT(K-l)*lE-9 
SUM=(0. 0,0. 0) 

DO 111 1=1, K 

SUM=SUM+Y4( I )*CONJG( Y4( I ) ) 

111 CONTINUE 

438 WRITE ( * j * ) 1 ***************************************** * 

WRITE(* **) ' ENTER RELATIVE SMOOTHING PARAMETER : ' 

READ (*,*) L9 
L1=CMPLX( L9 , 0. 0) 

L2=SUM/CMPLX(FLOAT(K+l) ,0. 0)*L1 
DO 1 1=1, K 

TUT( I )=Y5( I )*CONJG( Y4( I ) )*YSC( I ) /( Y4( I )*CONJG( Y4( I ) )+L2) 
1 CONTINUE 

TITLE=' SPECTRAL PLOT' 

YS=' DECONVOLUTION' 

XS=' X-AXIS FREQUENCY IN GHZ' 

NUMPT=K 

Q ************************************************* 

CALL MAGNITUD(TUT, NUMPT, MAG) 

Q * * * * * * * * * ***** * * * * * * * * * * * * * * * * ******************* 

XMI1=0. 0 

XMA1=1. O/T/2. 0*lE-9 

Q ************************************************* 

CALL P12(TITLE, NUMPT, XMI1,XMA1, MAG, XS, ID, YS,L9) 
c************************************************* 

CALL IFFF(K,T,TUT) 

Q************** *********************************** 

DO 222 1=1, K 
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R( I )=REAL( TUT( I ) ) 

222 CONTINUE 

XMI 1=0. 0 
XMA1=T2 

XS=' X-AXIS TIME IN NSEC' 
TITLE=' TIME SERIES PLOT' 
YS=' DECONVOLUTION' 



YN= ' N ' 

WRITE ( * j * ) 1 ***************************************************** 1 
WRITE(*’*) 'WOULD YOU LIKE TO SAVE THE DECONVOLUTION RESULT ?: 
READ(*,999) YN 

IF ( YN. EQ. 'Y'.OR. YN.EQ. ' y ' ) GO TO 4376 
GOTO 4388 



4376 WRITE ( ■* , * ) 1 ************************************** 1 



WRITE(*,*)' ENTER FILE NAME FOR THE TIME DOMAIN DATA: ' 

READ(*,999) NAME4 

CALL OUT( NAME4 , ID , TYPE , DATE , F 1 , SUB AV , T2 , DATAB , K , R ) 

4388 YN= ' Y ' 

WRITE ( * ^ * ) 1 *************************************************** 1 
WRITE (*’*) 'REPEAT WITH NEW SMOOTHING PARAMETER ? (Y/N): ' 

READ(* ,999) YN 

IF ( YN.EQ. 'Y'.OR. YN.EQ. 'y') GO TO 438 
999 FORMAT(A) 

RETURN 

END 





c ************************************************ 
C F = ARRAY OF Y VARIABLE 

C XS= $ X-SCALE (STRING) 

C YS= $ Y-SCALE (STRING) 

C ID= IDENTIFICATION OF THIS GRAPH 

C 



SUBROUTINE P 1 1( TITLE ,NPTS , XMINN , XMAXX , FFF , XS , ID , YS) 



c 



C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 



RM-FORTRAN Subroutine to Plot a Solid Line with TITLE 
Based on Program PLOT. FOR. Dec 1987 by M. A. Morgan 



INPUT DATA FORMAT 



TITLE - 64 Space Title 
NPTS - # Data Points 



XMIN - Real Min X value 
XMAX - Real Max X value 
F(N) - Input Data Array 
IHC - Hardcopy (l=Yes) 



CHARACTER*! YN,DUM 
CHARACTER* 6 4 TITLE, XS, ID, YS 
REAL X( 1025) ,F( 1025) , XMIN, XMAX, FFF( 1025) , XMINN, XMAXX 
INTEGER* 2 N , JROW , JCOL , I SYMBL , ITYPE , CYAN , WHITE , YELLOW 



IHC=0 
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WRITE (*,*) ' WOULD YOU LIKE TO HAVE A HARD COPY: ' 

READ(*, 100) YN 

IFCYN.EQ. 'Y'.OR.YN.EQ. V) IHC=1 

WHITE=7 

CYAN=11 

DO 1111 1=1 ,NPTS 
F( I )=FFF( I ) 

1111 CONTINUE 

XMIN=XMINN 

XMAX=XMAXX 

YELL0W=14 

N=NPTS 

DX=( XMAX-XMIN ) / ( NPTS - 1 . 0) 

FMIN=0. 0 
FMAX=0. 0 
DO 22 K=1,NPTS 
X(K)=XMIN+(K-1. 0)*DX 
IF(F(K). LT.FMIN) FMIN=F(K) 

IF(F(K). GT. FMAX) FMAX=F(K) 

22 CONTINUE 

IF(FMIN. GT. 0. 0) FMIN=0. 0 
IF(FMAX. LT. 0. 0) FMAX=0. 0 

C Computing Scale Factors for Vertical Axis 

ABSMIN=ABS(FMIN) 

ABSMAX=AB S ( FMAX ) 

YMAX=AMAX1( ABSMIN , ABSMAX) 

NSCL=INT(L0G10( YMAX ) ) 

IF (YMAX. LT. 1. 0) NSCL=NSCL-1 
YSCL=10. **NSCL 
FMIN=FMIN/YSCL 
FM AX=FM AX / YS C L 
ABSMIN=ABSMIN/YSCL 
ABSMAX=ABSMAX/YSCL 
DO 33 K=1 ,NPTS 

33 F(K)=F(K)/YSCL 

YMIN=0. 0 

IF(FMIN. EQ. 0. 0) GO TO 37 
35 YMIN=YMIN+0. 5 

IF( ABSMIN. GT. YMIN) GO TO 35 
YMIN=YMIN*FM IN/ ABSMIN 

37 CONTINUE 

YMAX=0. 0 

IF(FMAX. EQ. 0. 0) GO TO 41 
39 YMAX=YMAX+0. 5 

IF( ABSMAX. GT. YMAX) GO TO 39 
YMAX=YMAX*FMAX/ ABSMAX 

41 CONTINUE 

CALL QBEEP 

CALL QSMODE( 2) 

WRITE(*,*) 'Press RET to View Plot. .. RET Again to Clear Screen' 

IFCIHC.NE. 1) GO TO 42 

WRITE(*,*) 'A Hardcopy Will Be Made. .. Check that Printer is ON' 

42 CONTINUE 

READ(*,100) DUM 

C Calling GRAFMATIC Routines and Plotting Solid Line Graph 

ITYPE=1 
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ISYMBL=-2 
CALL QSMODE( 16) 

CALL QPTXT( 64, TITLE, YELLOW ,16,23) 

CALL QPTXT(64,ID,CYAN, 16,19) 

CALL QPTXT( 64, XS, CYAN, 16,20) 

CALL QPTXTC 64, YS, CYAN, 16,21) 

CALL QPLOT( 110, 540 , 20 , 240 , XMIN , XMAX , YMIN , YMAX , 0. ,0. ,0,1. ,1.5) 
CALL QSETUP(0, CYAN, ISYMBL, CYAN) 

XMAJOR=XMAX/5. 0 

CALL QXAXI S( XMIN, XMAX, XMAJOR, 1,1,2) 

YMAJ0R=1. 0 

CALL QYAXI S ( YMI N , YMAX , YMA JOR ,1,1,2) 
JR0W=14+(ABS(YMIN)/( ABS( YMAX) +ABS( YMIN) ) )*220 
JC0L=80-430*XMIN/( XMAX -XMIN) 

CALL QGTXT(3, '0. O' , WHITE, JC0L,JR0W,0) 

CALL QPTXTC 1, ‘s' , YELLOW, 5 , 16) 

CALL QPTXT( 1 , ' c ' , YELLOW ,5,15) 

CALL QPTXT(l,'a' , YELLOW, 5 , 14) 

CALL QPTXT( 1 , ' 1 ' , YELLOW ,5,13) 

CALL QPTXT( 1 , ' e ' , YELLOW ,5,12) 

CALL QPTXT( 1 , ' X ' YELLOW ,5,10) 

CALL QPTXT(2, ' 10 ' , YELLOW, 4, 8) 

CALL QPTXTC 1 , ' * ' , YELLOW ,5,7) 

CALL QPTXTC 1 , ' * ' , YELLOW ,5,6) 

CALL QCMOV( 0 ,6) 

WRITE (*, 150) NSCL 
CALL QTABL( ITYPE , N , X , F) 

READ(*,100) DUM 
IF(IHC.NE. 1) GO TO 44 
CALL PRTSC 

44 CONTINUE 

CALL QSM0DE(2) 

100 FORMAT(A) 

150 F0RMAT(4X,I3) 

RETURN 

END 







THIS SUBROUTINE PLOTS TWO GRAPHS ON THE SAME AXIS 



F = Y1 ARRAY VARIABLE ( BLUE) 

FI = Y2 ARRAY VARIABLE (GREEN) 

AFTER PLOT THIS PROGRAM BE AWARE OF THE F1,F 
IT WILL SENT BACK A SCALE F,F1 

XS = $ X-SCALE 
YS = $ Y-SCALE 
ID = IDENTIFICATION 



SUBROUTINE PL08( TITLE , NPTS , XMINN , XMAXX , PUM , XS , YS , PUM1 , ID) 



RM-FORTRAN Subroutine to Plot a Solid Line with TITLE 
Based on Program PLOT. FOR. Dec 1987 by M. A. Morgan 

INPUT DATA FORMAT 
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TITLE - 64 Space Title 
NPTS - # Data Points 
XMIN - Real Min X value 
XMAX - Real Max X value 
F(N) - Input Data Array 
IHC - Hardcopy (l=Yes) 

CHARACTER* 1 YN,DUM 
CHARACTER*64 TITLE ,XS , ID ,YS 

REAL X( 1025) ,F( 1025) ,F1( 1025) ,PUM( 1025) ,PUM1( 1025) ,XMINN,XMAXX 
I NTEGER*2 N , JROW , JCOL , I SYMBL, I TYPE , CYAN , WHITE , YELLOW , GREEN 

Q *** ************ ************** ******** ********** * 

WRITE (*,*) ’WOULD YOU LIKE TO HAVE A HARD COPY ?: ’ 

READ(*,100) YN 
IHC=0 

IF( YN. EQ. ' Y ’ . OR. YN. EQ. ’ y ' ) IHC=1 
DO 159 1=1, NPTS 
F(I)=PUM(I) 

Fl( I )=PUM1( I ) 

159 CONTINUE 

XMIN=XMINN 

XMAX=XMAXX 

(J******************'A ***************************** 

WHITE=7 

CYAN=11 

YELL0W=14 

GREEN=2 

N=NPTS 

DX=(XMAX-XMIN)/(NPTS-1. 0) 

FMIN=0. 0 
FMAX=0. 0 

FMIN1=0. 0 
FMAX1=0. 0 

DO 22 K=1 ,NPTS 
X(K)=XMIN+(K-1. 0)*DX 
IF(F(K). LT. FMIN) FMIN=F(K) 

IF(F(K) . GT. FMAX) FMAX=F(K) 

IF(F1(K). LT. FMIN1) FMIN1=F1(K) 

IF(F1(K). GT. FMAX1) FMAX1=F1(K) 

22 CONTINUE 

FMIN=AMIN1( FMIN, FMIN 1) 

FMAX=AMAX 1 ( FMAX , FMAX 1 ) 

IF(FMIN. GT. 0. 0) FMIN=0. 0 
IF(FMAX. LT. 0. 0) FMAX=0. 0 

C Computing Scale Factors for Vertical Axis 

ABSMIN=ABS(FMIN) 

ABSMAX=ABS ( FMAX ) 

YMAX=AMAX 1(ABSMIN, ABSMAX) 

NSCL=I NT( LOG 1 0 ( YM AX ) ) 

IF (YMAX. LT. 1. 0) NSCL=NSCL-1 
YSCL=10.**NSCL 
FMIN=FMIN/YSCL 
FMAX=FMAX/YSCL 
ABSMIN=ABSMIN/YSCL 
ABSMAX=ABSMAX/YSCL 
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DO 33 K=1 ,NPTS 

F1(K)=F1(K) /YSGL 

33 F(K)=F(K) /YSCL 

YMIN=0. 0 

IF(FMIN. EQ. 0. 0) GO TO 37 
35 YMIN=YMIN+0. 5 

IF(ABSMIN. GT. YMIN) GO TO 35 
YMIN=YMIN*FMIN/ABSMIN 

37 CONTINUE 

YMAX=0. 0 

IF(FMAX. EQ. 0. 0) GO TO 41 
39 YMAX=YMAX+0. 5 

IF(ABSMAX. GT. YMAX) GO TO 39 
YMAX=YMAX ,V FMAX/ABSMAX 

41 CONTINUE 

CALL QBEEP 

CALL QSMODE( 2) 

WRITE(*, ,V ) 'Press RET to View Plot. .. RET Again to Clear Screen' 

IF(IHC. NE. 1) GO TO 42 

WRITE( ,V ,*) 'A Hardcopy Will Be Made. .. Check that Printer is ON' 

42 CONTINUE 

READ(*, 100) DUM 

C Calling GRAFMATIC Routines and Plotting Solid Line Graph 

ITYPE=1 
ISYMBL=-2 
CALL QSM0DE( 16) 

CALL QPTXT( 64, TITLE, YELLOW, 16,23) 

CALL QPTXT(64,ID,GREEN, 16,19) 

CALL QPTXT( 64, XS, WHITE, 16,20) 

CALL QPTXT( 64, YS, CYAN, 16,21) 

CALL QPL0T(110, 540, 20, 240, XMIN,XMAX, YMIN, YMAX, 0. ,0. ,0,1. ,1.5) 
CALL QSETUP(0, CYAN, ISYMBL, CYAN) 

XMAJ0R=XMAX/5. 0 

CALL QXAXISf XMIN , XMAX , XMAJOR ,1,1,2) 

YMAJ0R=1. 0 

CALL QYAXI S ( YM I N , YMAX , YMA JOR ,1,1,2) 

JR0W=l4+( ABS( YMIN)/( ABS( YMAX)+ABS( YMIN) ) ) ,v 220 
JC0L=80-430*XMIN/ (XMAX-XMIN) 

CALL QGTXT(3, '0. O' , WHITE, JC0L,JR0W,0) 

CALL QPTXT(1, 's' , YELLOW, 5 , 16) 

CALL QPTXT( 1 , ' c ' , YELLOW ,5,15) 

CALL QPTXT( 1 , ' a ' , YELLOW ,5,14) 

CALL QPTXT( 1, ' l' , YELLOW, 5, 13) 

CALL QPTXT(l,'e' , YELLOW, 5 , 12) 

CALL QPTXT(1,'X' , YELLOW, 5 , 10) 

CALL QPTXT( 2 , ' 10 ' , YELLOW ,4,8) 

CALL QPTXT( 1 , ' * ' , YELLOW ,5,7) 

CALL QPTXT(1, '*' , YELLOW, 5, 6) 

CALL QCMOV( 0 ,6) 

WRITE(*, 150) NSCL 
CALL QTABL( ITYPE ,N,X,F) 

QJctcfcfrjc'kMcjcJc'k'tcicfcrk'kMcicjtjcMe&jcjtJr&itMeicic'kMcWc 



C PLOT ANOTHER GRAPH 
GREEN=10 



ISYMBL=-2 
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000000000 0000 ooooo 



ITYPE=1 

CALL QSETUP(0, GREEN, ISYMBL, GREEN) 
CALL QTABL( ITYPE,N,X,F1) 
READ(*, 100) DUM 
IF( IHC. NE. 1) GO TO 44 
CALL PRTSC 

44 CONTINUE 

CALL QSMODE( 2) 

100 FORMAT(A) 

150 FORMAT( 4X, 13) 

RETURN 

END 



Qic'kjcjrtrtrfrtelcjcjcjcjn'cft ********Vf9V*VoV>V***^r********^V****** 

F = ARRAY OF Y VARIABLE 
XS= $ X-SCALE (STRING) 

YS= $ Y- SCALE (STRING) 

ID= IDENTIFICATION OF THIS GRAPH 

SUBROUTINE PI 2( TITLE , NPTS , XM INN , XMAXX , FFF , XS , ID , YS , L9 ) 

RM-FORTRAN Subroutine to Plot a Solid Line with TITLE 



Based on Program PLOT. FOR. Dec 


1987 by M. A. Morgan 


INPUT DATA FORMAT 




TITLE 


- 64 Space Title 


NPTS 


- # Data Points 


XMIN 


- Real Min X value 


XMAX 


- Real Max X value 


F(N) 


- Input Data Array 


IHC 


- Hardcopy (l=Yes) 



CHARACTER* 1 YN , DUM 
CHARACTER*64 TITLE , XS , ID , YS , RELAT 
REAL X( 1025) ,F( 1025) ,XMIN,XMAX,FFF( 1025) ,XMINN, XMAXX, L9 
INTEGER* 2 N , JROW , JCOL , ISYMBL, I TYPE , CYAN , WHITE , YELLOW 

IHC=0 

WRITE(*,*) ' WOULD YOU LIKE TO HAVE A HARD COPY: ' 

READ(*, 100) YN 

IF(YN. EQ. 'Y'.OR.YN.EQ. V) IHC=1 
WHITE=7 
CYAN=11 

DO 1111 1=1, NPTS 

F( I)=FFF( T) 

1111 CONTINUE 

XMIN=XMINN 

XMAX=XMAXX 

YELL0W=14 

N=NPTS 

DX=(XMAX-XMIN)/(NPTS-1. 0) 

FMIN=0. 0 
FMAX=0. 0 
DO 22 K=1 ,NPTS 
X(K)=XMIN+(K-1. 0)*DX 
IF(F(K). LT. FMIN) FMIN=F(K) 
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IF(F(K). GT. FMAX) FMAX=F(K) 

CONTINUE 

IF(FMIN. GT. 0. 0) FMIN=0. 0 
IF(FMAX. LT. 0. 0) FMAX=0. 0 
Computing Scale Factors for Vertical Axis 
ABSMIN=ABS(FMIN) 

ABSMAX=ABS(FMAX) 

YMAX=AMAX1( ABSMINjABSMAX) 

NSCL=INT( LOG 1 0 ( YMAX ) ) 

IF (YMAX. LT. 1. 0) NSCL=NSCL-1 
YSCL=10. **NSCL 
FMIN=FMIN/YSCL 
FM AX=FM AX / Y S C L 
ABSMIN=ABSMIN/YSCL 
ABSMAX=ABSMAX/YSCL 
DO 33 K=1 ,NPTS 
F(K)=F(K)/YSCL 
YMIN=0. 0 

IF(FMIN.EQ. 0. 0) GO TO 37 
YMIN=YMIN+0. 5 

IF(ABSMIN.GT.YMIN) GO TO 35 
YMIN=YMIN*FMIN/ABSMIN 
CONTINUE 
YMAX=0. 0 

IF(FMAX.EQ. 0. 0) GO TO 41 
YMAX=YMAX+0. 5 

IF(ABSMAX. GT. YMAX) GO TO 39 
YMAX=YMAX*FMAX/ ABSMAX 
CONTINUE 
CALL QBEEP 
CALL QSMODE( 2) 

WRITE(*,*) 'Press RET to View Plot. ..RET Again to Clear Screen 

IF( IHC. NE. 1) GO TO 42 

WRITE(*,*) 'A Hardcopy Will Be Made. .. Check that Printer is ON 

CONTINUE 
READ(*,100) DUM 

Calling GRAFMATIC Routines and Plotting Solid Line Graph 

ITYPE=1 
ISYMBL=-2 
CALL QSM0DE( 16) 

CALL QPTXT( 64 , TITLE , YELLOW ,16,23) 

CALL QPTXT( 64 , ID , CYAN ,16,19) 

CALL QPTXT( 64, XS, CYAN, 16,20) 

CALL QPTXT( 64, YS, CYAN, 16,21) 

RELAT=' RELATIVE SMOOTH =' 

CALL QPTXT(64,RELAT,CYAN, 16,18) 

CALL QQNPUT( 265 , 253 , L9 , 3) 

CALL QPL0T( 1 10 , 540 , 20 , 240 , XMIN , XMAX , YMIN , YMAX , 0. ,0. ,0,1. ,1.5) 
CALL QSETUP( 0, CYAN, ISYMBL, CYAN) 

XMAJOR=XMAX/5. 0 

CALL QXAXIS ( XMIN , XMAX , XMA JOR ,1,1,2) 

YMAJ0R=1. 0 

CALL QYAXI S( YM IN , YMAX , YMA JOR ,1,1,2) 
JR0W=14+(ABS(YMIN)/(ABS(YMAX)+ABS(YMIN)))*220 
JC0L=80 -430*XMIN/ ( XMAX -XMIN ) 

CALL QGTXT(3, '0. O' , WHITE, JCOL,JROW,0) 
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CALL QPTXT( 1, ' s' , YELLOW, 5 , 16) 
CALL QPTXT( 1 , ' c ' , YELLOW ,5,15) 
CALL QPTXT( 1, ' a' , YELLOW, 5 , 14) 
CALL QPTXT( 1 , * 1 ' , YELLOW, 5 ,13) 
CALL QPTXT ( 1 , ' e ' , YELLOW ,5,12) 
CALL QPTXT(1,'X' YELLOW, 5, 10) 
CALL QPTXT ( 2 , ' 10 f , YELLOW ,4,8) 
CALL QPTXT( 1 , ' * ' , YELLOW ,5,7) 
CALL QPTXT(1, , YELLOW, 5, 6) 
CALL QCMOV( 0,6) 
WRITE(*, 150) NSCL 
CALL QTABL( ITYPE,N,X,F) 
READ(*, 100) DUM 
IF( IHC. NE. 1) GO TO 44 
CALL PRTSC 

44 CONTINUE 

CALL QSM0DE(2) 

100 FORMAT( A) 

150 F0RMAT(4X,I3) 

RETURN 

END 



SUBROUTINE PLO 7 7 ( Y 1 , Y5 , T2 , NUMPT , ID 1 ) 

0*** *** ******** ** ** ***** * * *** ***** * * *** ** * ** * *** ** 

Y2 = Y5 = BACKGROUND DATA ARRAY * 

Y1 = SCATTERED DATA ARRAY * 

0 *** ***** *** *** *** ****** **** *** ** **** ***** ******** 

INTEGER NUMPT , DATAB , SUBAV , FI , NSHIF , SET 

REAL Yl( 1024) , Y2( 1024) , Y3( 1024) ,Y4( 1024) , Y5( 1024) 

REAL XMIN,XMAX 

CHARACTER*^ NAME 1 , NAME 2 , NAME 3 , QN 
CHARACTER* 1 YN 

CHARACTER*64 YS , XS , ID , TITLE , DATE 
CHARACTER*64 ID 1 , TYPE 1 
CHARACTER* 6 4 ID2,TYPE2 
CHARACTER*64 ID3,TYPE3 
REAL TIME ,T2 

WR I TE ( * * ) f **************************************************** 1 
WRITE(*,’*) ’SHIFTING SCATTERED WAVEFORM & SUBTRACTING BACKGROUND' 
WR I TE ( " ) ’ ** yr**********^'*********** ********* *** Vf** ***** 1 

XMIN=0. 6 



XMAX=T2 

DO 22 K=l, NUMPT 
22 Y2(K)=Y5(K) 

SET=0 

Q*****************Vf****************************** 



33 CALL SHIFT(TITLE , YS ,XS , ID) 

Q**************'******************’**************** 



CALL PL08 ( TITLE , NUMPT , XMIN , XMAX , Y1 , XS , YS , Y2 , ID ) 

Q********************************Vc *************** 

C SHIFT THE CURVE 

Q************************************************ 



CALL NSHIFT( NSHIF) 

IF( NSHIF. EQ. 0) GO TO 40 

CALL G00D(Y1,Y2, NSHIF, NUMPT, SET) 

GO TO 33 
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40 CONTINUE 

Q************************************************ 

C SUBTRACT BACKGROUND FROM THE WAVEFORM **** 

DO 44 1=1 ,NUMPT 
44 Y3(I)=Y1(I)-Y2(I) 

CALL SUBTR( TITLE, XS,YS) 

CALL P 11( TITLE ,NUMPT,XMIN,XMAX, Y3 ,XS , ID1 , YS) 

WRITE(*, ,V ) 'DO YOU WANT TO RE-SHIFT AND SUBTRACT ? (Y/N): 
READ(* , 300) YN 

IF(YN. EQ. 'y'.OR.YN. EQ. ' Y* ) GO TO 33 
DO 77 1=1 ,NUMPT 
77 Yl( I )=Y3( I ) 

300 FORMAT(A) 

RETURN 

END 

Q ********* ** *** ********* * ************************* 

SUBROUTINE NSHIFT(NSHIF) 

INTEGER NSHIF 

WRITE ( * , ,f ) 1 ****************************************** 1 
WRITEOVO 

WRITE(*,*) 'SCATTERED WAVEFORM SHIFTING: ' 

WRITE(*,*) 

WRITE(*,*) ' +n shifts scat waveform to the RIGHT n-points' 
WRITE(*,*) ' -n shifts scat waveform to the LEFT n-points' 
WRITE(*,*) ' 0 DISPLAYS subtracted waveform ' 

WRITEOVO 

WRITEC*,*) 'ENTER SHIFT (+/-INTEGER) : ' 

READ(*,*) NSHIF 

RETURN 

END 



C ***** 



:**********************■************** 



SUBROUTINE SUBTR( TITLE ,XS ,YS) 
CHARACTER*64 TITLE ,XS , ID, YS 
TITLE=' SUBTRACTED WAVEFORM' 
XS=' X-AXIS TIME IN NS' 

YS=' Y-AXIS MAGNITUDE IN VOLT' 

RETURN 

END 



Q * * * * * * * *5 V * * * * * * * * * 7 V * * * * * * * * * * * * * * * * * * * * * * ******** 



SUBROUTINE SHIFT(TITLE,YS,XS,ID) 
CHARACTER*64 TITLE ,YS ,XS , ID 
TITLE=' TIME SERIES PLOT 1 
YS=' BLUE -TARGET' 

I D= ' GREEN - BACKGROUND ' 

XS=' SHIFT-TARGET ONLY' 

RETURN 

END 



c** ***** *** ******** * 



***************************** 
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SUBROUTINE GOOD( Y1 , Y2 ,NSHIF,NUMPT, SET) 
PROGRAM TO SHIFT Y2 CURVE 



Y1 TO THE LEFT 
Y1 TO THE RIGHT 

IS THE SHIFTING CURVE 
IS THE BACKGROUND CURVE 
IS THE NUMBER OF POINT TO SHIFT 
IS THE SUM OF THE POINT TO SHIFT 
IN ORDER TO ZERO PADDING 
IS THE NUMBER OF POINT 



* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 



C 
C 

c 

C - SHIFT 

C + SHIFT 

C 

C Y1 

C Y2 

C NSHIF 

C SET 

C 

C NUMPT 

C 

REAL Yl( 1024), Y2( 1024), Y3( 1024) 

INTEGER NSHIF, NUMPT, SET 

5 IF( NSHIF. LT. 0 ) THEN 

NSHIF=-1*( NSHIF) 

DO 20 I=NSHIFT,NUMPT-NSHIF-SET 
Y3( I ) = Y 1 ( I+NSHIF) 

20 CONTINUE 

DO 10 I=NUMPT-NSHIF-SET, NUMPT 
Y2( I)=0. 0 
Y3( I)=0. 0 

10 CONTINUE 

DO 15 1=1 ,NSHIF+SET+1 
Y2( I)=0. 0 
Y3( I )=0. 0 

15 CONTINUE 

ELSE IF( NSHIF. GT. 0 ) THEN 
N=NSHIF+SET+1 
DO 40 I=N, NUMPT 

Y3( I)=Y1( I-NSHIF) 

40 CONTINUE 

DO 30 1=1, N 

Y2( I )=0. 0 
Y3( I )=0. 0 

30 CONTINUE 

DO 35 I=NUMPT-NSHIF-SET, NUMPT 
Y2( I )=0. 0 
Y3( I)=0. 0 

35 CONTINUE 

ELSE 
END IF 

SET=SET+NSHIF 
DO 50 1=1, NUMPT 
Yl( I)=Y3( I) 

50 CONTINUE 
RETURN 
END 



C 

c 



SUBROUTINE OUT( NAME , ID ,TYPE , DATE ,F 1 , SUBAV , T2 , DATAB , NUMPT , Y ) 
ID $ = CHARACTER $ TARGET IDENTIFICATION 
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oooooooo 



TYPE $ = TYPE OF WAVEFORM 

DATE $ = DATE OF EXPERIMENT TAKE PLACE 

FI = DUMMY STRING 

SUBAV = SUB AVERAGE 

NUMPT = NUMPT OF POINT 

DATAB = DATA BLOCK 

Y = Y ARRAY 

T2 = TIME WINDOW 

CHARACTER*64 ID , TYPE , DATE ,F1 

CHARACTER* 16 NAME 

REAL Y( 1024) ,T2 

INTEGER SUBAV, NUMPT, DATAB 

OPEN( UNIT* 10 ,FILE=NAME , STATUS* ' UNKNOWN ’ ) 

WRITE (10,*) ID 

WRITE(10,*) TYPE 

WRITE (10,*) DATE 

WRITEUO,*) FI 

WRITE( 10,22) SUBAV 

WRITE( 10,21) T2 

WRITE( 10,22) DATAB 

WRITE( 10,22) NUMPT 

DO 10 1=1, NUMPT 

WRITE( 10,20) Y( I) 

10 CONTINUE 

20 FORMAT(E12. 6) 

21 FORMAT(F5. 2) 

22 FORMAT( 15) 

30 FORMAT(A) 

CLOSE(IO) 

RETURN 

END 
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APPENDIX B. THE ACQUISITION QUICK BASIC PROGRAM 



The following computer program acquires the measurement data from the HP 
54120T as described in chapter 5. The program is written in QUICK BASIC 



1000 

1020 

1040 

1060 

1080 

1240 

1321 

1322 

1324 

1325 

1326 

1327 

1328 

1329 

1330 

1341 

1342 

1343 

1344 

1345 

1346 

1347 
1360 
1380 
1400 
1420 
1440 
1460 
1480 
1500 
1520 
1540 
1560 
1580 
1600 
1620 
1640 
1660 
1680 
1700 
1720 
1740 
1760 
1780 
1800 
1820 
1840 



REM ******************************************* 
REM * MODIFY 4/22/88 * 

WRITTEN BY LT. JG SOONPUEN SOMAPEE * 

ACQUISITION * 

* 



REM * 

REM * 

REM * 

REM * * 

REM ******************************************* 

D$=DATE$ 

DUMMY$=TIME$: PRINT ’’DATE ";D$: PRINT "TIME DUMMY? 

COLOR 7,1 

PRINT "WELCOME TO ACQUISITION PROGRAM" 

PRINT "THERE ARE 3 TYPE OF MEASUREMENT" 

PRINT "1. BACKGROUND WAVETYPE" 

PRINT "2. CALIBRATE WAVETYPE" 

PRINT "3. TARGET WAVETYPE" 

PRINT " PRESS THE FOLLOWING KEY " 

PRINT " L ( USED DATA LAST RESPONSE) " 

PRINT " D ( SYSTEM DEFAULT ) " 

V$=INKEY$: IF V$="" THEN 1344 
IF V$="L" OR V$ ="1" THEN 1960 
IF V$="d" OR V$ ="D" THEN 1360 
GOTO 1344 

REM **** *** ** ******* ********** ******* * ********* 

REM * DEFAULT * 

REM ** ***************************************** 

OPEN "DATA.DAT" FOR OUTPUT AS #1 
Pl=512 
PRINT #1,P1 
N=1024 
PRINT# 1,N 
ND=1 

PRINT# 1 ,ND 

T$="4 inch Thin wire" 

PRINT# 1,T$ 

D$=DATE$ 

PRINT# 1,D$ 

T2=l. 563E-10 
PRINT# 1,T2 

W$="Background Wave Type" 

PRINT# 1,W$ 

Fl=20 

PRINT# 1, FI 
Sl=16 

PRINT# 1, SI 
G1=S1 

DSE$="Backg" 

PRINT# 1,DSE$ 
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1860 DUMMY$=TIME$ 
1880 PRINT# 1, DUMMY? 
1900 DELAY1=38 
1920 PR INT# 1 , DELAY 1 
1940 CLOSE# 1 
1960 RE M 



1980 


REM Optionally include the following declarations 




2000 


REM in your program. 




2020 


REM They provide 




appropriate mnemonics by which 




2040 


REM to reference 




commonly used values. 




2060 


REM Some mnemonics (GET%, ERR%, 




2080 


REM END%, ATN%) 


are preceded by B in order to 




2100 


REM distinguish them from 




2120 


REM BASICA keywords. 




2140 


REM 








2160 


REM GPIB Commands 






2180 


UNL% = &H3F 


1 


GPIB unlisten command 




2200 


UNT% = &H5F 


1 


GPIB untalk command 




2220 


GTL% = &H1 


1 


GPIB go to local 




2240 


SDC% = &H4 


1 


GPIB selected device clear 




2260 


PPC% = &H5 


1 


GPIB parallel poll configure 




2280 


BGET% = &H8 


1 


GPIB group execute trigger 




2300 


TCT% = &H9 


1 


GPIB take control 




2320 


LL0% = &H11 


1 


GPIB local lock out 




2340 


DCL% = &H14 


1 


GPIB device clear 




2360 


PPU% = &H15 


1 


GPIB ppoll unconfigure 




2380 


SPE% = &H18 


t 


GPIB serial poll enable 




2400 


SPD% = &H19 


f 


GPIB serial poll disable 




2420 


PPE% = &H60 


1 


GPIB parallel poll enable 




2440 


PPD% = &H70 


t 


GPIB parallel poll disable 




2460 


REM 








2480 


REM GPIB status 


bit vector 




2500 


REM global variable IBSTA% and wait mask 




2520 


BERR% = &H8000 


1 


Error detected 




2540 


TIM0% = &H4000 


1 


Timeout 




2560 


BEND% = &H2000 


1 


E0I or EOS detected 




2580 


SRQI% = &H1000 


t 


SRQ detected by CIC 




2600 


RQS% = &HC00 


! 


Device needs service 




2620 


CMPL% = &H100 


1 


1/0 completed 




2640 


L0K% = &H80 


1 


Local lockout state 




2660 


REM% = 8H40 


1 


Remote state 




2680 


CIC% = &H20 


1 


Cont roller- In -Charge 




2700 


BATN% = &H10 


1 


Attention asserted 




2720 


TACS% = 6cH8 


1 


Talker active 




2740 


LACS% = &H4 


t 


Listener active 




2760 


DTAS% = &H2 


1 


Device trigger state 




2780 


DCAS% = &H1 


1 


Device clear state 




2800 


REM 








2820 


REM Error messages returned in global variable IBERR% 


2840 


EDVR% = 0 


1 


DOS error 




2860 


ECIC% ^ 1 


T 


Function requires GPIB-PC to be 


CIC 


2880 


EN0L% = 2 


1 


Write function detected no Listeners 


2900 


EADR% = 3 


t 


Interface board not addressed correctly 


2920 


EARG% = 4 


t 


Invalid argument to function call 


2940 


ESAC% = 5 


1 


Function requires GPIB-PC to be 


SAC 


2960 


EAB0% = 6 


1 


I/O operation aborted 
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2980 


ENEB% = 7 


' Non-existent interface board 


3000 


EOIP% = 10 


1 1/0 operation started before 


3020 


REM 


operation completed 


3040 


ECAP% = 11 


No capat 


•ility for operation 


3060 


EFSO% = 12 


File system operation error 


3080 


EBUS% = 14 


Command 


error during device i 


3100 


ESTB% = 15 


Serial poll status byte lost 


3120 


ESRQ% =16 


' SRQ remains asserted 


3140 


REM 






3160 


REM EOS mode bits 




3180 


BIN% = &H1000 


Eight bit compare 


3200 


XE0S% = &H800 


' Send E0I 


! with EOS byte 


3220 


REOS% = &H400 


1 Terminate read on EOS 


3240 


REM 






3260 


REM Timeout values and meanings 


3280 


TN0NE% = 0 


Infinite 


s timeout (disabled) 


3300 


T10US% = 1 


' Timeout 


of 10 us (ideal) 


3320 


T30US% = 2 


Timeout 


of 30 us (ideal) 


3340 


T100US% = 3 


' Timeout 


of 100 us (ideal) 


3360 


T300US% = 4 


' Timeout 


of 300 us (ideal) 


3380 


T1MS% = 5 


Timeout 


of 1 ms (ideal) 


3400 


T3MS% = 6 


' Timeout 


of 3 ms (ideal) 


3420 


T10MS% = 7 


T imeout 


of 10 ms (ideal) 


3440 


T30MS% = 8 


Timeout 


of 30 ms (ideal) 


3460 


T100MS% = 9 


Timeout 


of 100 ms (ideal) 


3480 


T300MS% = 10 


Timeout 


of 300 ms (ideal) 


3500 


T1S% = 11 


Timeout 


of 1 s (ideal) 


3520 


T3S% = 12 


Timeout 


of 3 s (ideal) 


3540 


T10S% = 13 


Timeout 


of 10 s (ideal) 


3560 


T30S% = 14 


Timeout 


of 30 s (ideal) 


3580 


T100S% = 15 


Timeout 


of 100 s (ideal) 


3600 


T300S% = 16 


Timeout 


of 300 s (ideal) 


3620 


T1000S% = 17 


' Timeout 


of 1000 s (maximum) 


3640 


REM 






3660 


REM Miscellaneous 




3680 


S% = &H8 


' Parallel 


. Poll sense bit 


3700 


LF% = &HA 


1 Line feed character 


3720 


REM 






3740 


REM Application 


program variables passed to 


3760 


REM GPIB functions 




3780 


REM 






3800 


CMD$ = SPACE$( 10) ' 


command buffer 


3820 


RD$ = SPACE$( 255) 


read data buffer 


3840 


WRT$ = SPACE$( 255) ' 


write data buffer 


3860 


BNAME$ = SPACE $( 7) ' 


board name buffer 


3880 


BDNAME$ = SPACE$( 7 ) ' 


board or device name 


3900 


FLNAME$ = SPACE$( 50) ' 


file name buffer 



3940 CLS 

3960 DIM W( 1500), RRR( 1500) 

3980 DIM RDSS$( 1024) 

4000 SCREEN 0 
4020 COLOR 7,1,3 

4040 OPEN "C: DATA. DATA" FOR INPUT AS #1 
4060 BEEP: CLS 

4080 REM **************************************** 
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4100 

4120 

4140 

4160 

4180 

4200 

4220 

4240 

4260 

4280 

4300 

4320 

4340 

4360 

4380 

4400 

4420 

4440 

4500 

4520 

4540 

4560 

4580 

4600 

4620 

4640 

4660 

4680 

4700 

4720 

4740 

4760 

4780 

4800 

4820 

4840 

4860 

4880 

4900 

4920 

4940 

4960 

4980 

5000 

5020 

5040 

5060 

5080 

5100 

5120 

5140 

5160 

5180 

5200 

5220 

5240 



REM * BEGIN MAIN PROGRAM * 

REM 



LOCATE 5,8: PRINT" DATA ACQUISITION 



REM 


************************************************************ 


REM 


* 


T$ 


= 


TARGET IDENTIFICATION 


* 


REM 


* 


w$ 


= 


WAVEFORM TYPE 


* 


REM 


* 


D$ 


= 


DATE 


* 


REM 


* 


N 


= 


NUMBER OF SUB -AVERAGE / WAVEFORM OR COUNT 


* 


REM 


* 


T2 


= 


DPO SAMPLING INTERVAL (nsec) 


* 


REM 


* 


SI 


= 


DPO MAXIMUM DPO VERTICAL SCALE(M) 


* 


REM 


* 


ND 




NUMBER OF DATA BLOCKS 


* 


REM 


* 


PI 


= 


NUMBER OF POINT / WAVEFORM 


* 


REM 


* 


Fl 


= 


DPO TIME WINDOW 


* 


REM 


* 


DSE$ 


= 


FILE NAME OF THE DATA OUTPUT FILE 


* 


REM 


* 


DEV$ 


= 


NAME OF THE HPDEV IN THE CONFIGURATION FILE* 


REM 


* 


DSO% 


zz 


THE NUMBER ASSOCIATED WITH THE DEVICE 


* 


REM 


* 


VEW$ 


zz 


STRING$ VARIABLE FOR SENDING TO HP54120 


* 


REM 


* 


DUMMY $ 


= 


TIME SET 


* 


REM 

REM 


************************************************************ 



INPUT#1 ,P1 

PRINT "NUMBER OF POINT "j PI 

INPUT#1 ,N 

PRINT "NUMBER OF SUB -AVERAGE( count) "jN 

INPUT#1,ND 

PRINT "NUMBER OF DATA BLOCK ND 

INPUT# 1 ,T$ 

PRINT "TARGET ID ";T$ 

INPUT#1,D$ 

PRINT "DATE ";D$ 

INPUT# 1 , T2 

PRINT "DPO SAMPLING INTERVAL ";T2 

INPUT#1 , W$ 

PRINT "WAVEFORM TYPE ";W$ 

INPUT# 1 ,F1 

PRINT "DPO TIMEWINDOWC NSEC ) FI 

INPUT# 1 , SI 

PRINT "DPO VERTICAL SCALE (M) ";S1 

G1=S1 



INPUT# 1, DSE$ 

PRINT "DATA FILE NAME OUTPUT ";DSE$ 

I NPUT# 1 , DUMMY $ 

PRINT "AUTO TIME SET DUMMY $ 

I NPUT# 1, DELAY 1 

PRINT "DELAY ( NSEC ) DELAY 1 

CLOSE# 1 

PRINT " 



COLOR 7,1: PRINT " PRESS THE FOLLOWING KEY 

PRINT " Y ( CHANGE THE DATA. DAT ) 

PRINT " N (NO CHANGE ) 

V$=INKEY$: IF V$="" THEN 5140 
IF V$="Y" OR V$ ="y" THEN 11260 
IF V$="N" OR V$ ="n" THEN 5220 



GOTO 5100 

dev$="hpdev" 

RE M ********** 



■** ******** 



************* ******* * 



fl 
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5260 REM * OPEN DEVICE AND RETURN THE UNIT * 

5280 REM * DESCRIPTOR ASSOCIATED WITH THE GIVEN * 
5300 REM * NAME * 

5320 REM ***************************************** 
5340 CALL IBFIND(DEV$ ,DSO%) 

5360 IF DSO% < 0 THEN PRINT "ERROR IN IBFIND" 

5380 RE M ***********************************■)>***** 
5400 REM * CLEAR SPECIFIED DEVICE * 

5420 REM ***************************************** 
5440 CALL IBCLR(DS0%) 

5460 REM ***************************************** 
5480 GOSUB 6080 
5500 CLS: BEEP 

5520 PRINT W$+" DATA STORAGE COMPLETE" 

5540 PRINT " MEASURE OTHER WAVEFORM TYPE(Y/N) " 

5560 Q$=INKEY$: IF Q$=" " THEN GOTO 5560 

5580 IF Q$="N" OR Q$="n" THEN GOTO 5640 

5600 IF Q$="Y" OR Q$="y" THEN GOTO 4040 

5620 GOTO 5560 

5640 STOP 

5660 RE M****************************************** 
5680 REM* SUBPROGRAM * 

5700 REM* * 

5720 REM * * **** **** ** **** **** ** ************ ** *** * * 



5740 


REM 


* 


TEXT REFERENCE 


* 


5760 


REM 


* 


"GURU II GPIB USER'S RESOURCE" 


* 


5780 


REM 


* 


"UTILITY FOR IBM PC" 


Vf 


5800 


REM 


* 


"TEKTRONIX 1986" 


* 



5820 REM ***************************************** 
5840 REM ***************************************** 
5860 REM * WRITE DATA FROM STRING SUBROUTINE * 
5880 REM ***************************************** 
5900 CALL I BWRT ( DSO% , VEW $ ) 

5920 RETURN 

5940 REM ***************************************** 



5960 


REM 


•J- 

4\ 


TEXT REFERENCE 


* 


5980 


REM 


* 


"GETTING STARTED GUIDE" 


* 


6000 


REM 


* 


"HP 54120T DIGITIZING " 




6020 


REM 


* 


"OSCILLOSCOPE" 


* 


6040 


REM 


* 


"HP COMPANY 1987" 


* 



6060 REM **** ************ *** ********************* * 
6080 REM ** *************************************** 
6100 REM * SETUP THE DISPLAY * 

6120 REM ***************************************** 
6140 VEW$="*RST" 

6160 GOSUB 5840 

6 180 REM ***************************************** 
6200 REM * CRT FORMAT 1 PROVIDE ONE DISPLAY AREA * 
6220 REM * AND USE 8 DIVISIONS FOR THE FULL * 

6240 REM * SCALE RANGE * 

6260 REM ***************************************** 
6280 VEW$=": DISPLAY: GRATICULE FRAME; FORMAT l" 

6300 GOSUB 5840 
6320 REM * 

6340 REM * TURN TDR STEP GENERATOR ON * 

6360 REM icicJrft'fc Vr *V V: * V? Vr V? VcVc Vc Vc Vr VcVcVr ic VoWr * ’sWcVc* Vc Vc* * Vc* VcVc 
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6380 

6400 

6420 

6440 

6460 

6480 

6500 

6520 

6540 

6560 

6580 

6600 

6620 

6640 

6660 

6680 

6700 

6720 

6740 

6760 

6780 

6840 

6860 

6880 

7020 

7040 

7060 

7080 

7100 

7120 

7140 

7160 

7180 

7200 

7220 

7240 

7260 

7280 

7300 

7320 

7340 

7360 

7380 

7400 

7420 

7440 

7460 

7480 

7500 

7520 

7540 

7560 

7580 

7600 

7620 

7640 



VEW$=": NETWORK: REFLECTION: PRESET" 

GOSUB 5840 

REM * *** ** * * * ** * ** * * * ** * * ****** ***** ****** * * * 

REM * PAGE 6-13 TURN ON CHANNEL 4 * 

REM ** **** * * * * * *** * *** **** * * * * *** ********* * 

VEW$=": VIEW CHANNEL4" 

GOSUB 5840 

REM ***************************************** 

REM * PAGE 6-4 TURN OFF CHANNEL1 * 

REM *** * * **** ******* ** ** * * * * * **** ** ** *** * * *** 

VEW$=": BLANK CHANNELl" 

GOSUB 5840 

REM ****** **** ***** ************** * ** ********* 

REM * SET UP CHANNEL RANGE AND OFFSET * 

REM * PAGE 17-6 * 

REM * TO BE 20 NANOSECONDS * 

REM * ** ***** ***** ******** *** ** ********** *** ** 

VEW$=": TIMEBASE: RANGE 20 NS" 

GOSUB 5840 

REM PR I NT 1 T ************************************************** 1 1 
REM PRINT "SET THE VERTICAL SCALE RANGE FROM 8 TO 640 MV 
VERT=S 1 

VEW$=": CHANNEL4: RANGE "+STR$(VERT)+"M" 

GOSUB 5840 

REM ***************************************** 

REM * TURN ON THE GRATITUDE * 

REM ******** *** ** *** * ** * ***** ** ** * ****** ***** 

VEW$=": DISPLAY: GRATICULE GRID" 

GOSUB 5840 

VEW$=": ACQ: BAND LOW" 

GOSUB 5840 

REM ** * * ** * ** * ** ** * * ** ** *** * ** * * ** ***** ****** 

REM * AVERAGING ON * 

T) I ’ V/ «.*» mJ «kV» J« fc-V «JS* •V* »*■» «|V *1* m)* 

k H M *% «« * w* *\ «> «> «« *\ /t <» /V m m A o o *\ A A 0 \ A n 7»7i A a a#v7w»7v7va a a 

VEW$=": ACQUIRE: TYPE AVERAGE" 

GOSUB 5840 

REM ***************************************** 

REM * ASKING FOR BANDWIDTH * 

REM 

p J ^ I ^ ^ ^ Jj. ^ ^ ^ ^ ^ ^ * 

PRINT "THE DEFAULT BANDWIDTH IS 12. 4GHZ " 

PRINT "WOULD YOU LIKE TO CHANGE TO 20 GHZ (Y/N)" 

Q$=INKEY$: IF Q$=" " THEN GOTO 7380 
IF Q$="N" OR Q$="n" THEN GOTO 7500 
IF Q$="Y" OR Q$="y" THEN GOTO 7460 
GOTO 7380 

VEW$=": ACQ: BAND HIGH" 

GOSUB 5840 
RD$=SPACE$( 255 ) 

VEW$=" ACQUIRE: COUNT "+STR$(N) 

GOSUB 5840 

VEW$=": ACQUIRE: POINTS "+STR$(P1) 

GOSUB 5840 

REM ********** **** ******* ***************** ** 

REM * SET TIME DELAY * 

REM * 
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7660 

7680 

7700 

7720 

7740 

7760 

7780 

7800 

7820 

7840 

7860 

7880 

7900 

7920 

7940 

7960 

7980 

8000 

8020 

8040 

8320 

8340 

8360 

8380 

8400 

8420 

8560 

8600 

8620 

8640 

8660 

8680 

8700 

8720 

8740 

8760 

8780 

8800 

8820 

8840 

8860 

8880 

8900 

8920 

8940 

8960 

8980 

9000 

9020 

9040 

9060 

9080 

9100 

9120 

9140 

9160 



REM **************************************** 

REM 

REM 

REM 

DELAY=DELAY1 
DUM?=STR?( DELAY) 

VEW?=": TIM: DEL f '+DUM$+"NS" 

GOSUB 5840 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM **************************************** 

REM HP 54120T SENT THE LAST BYTE RESPONSE "LF" ASH DECIMAL 10 
VEW$=": SYSTEM: HEADER OFF; : EOI ON" 

GOSUB 5840 

PRINT " * *************** ** *** *** ****** *** * ** *** *** *** ** * ** U 
REM **************************************** 

PRINT "I AM GOING TO DIGITIZE CHANNEL(Y/N)" 

Q?=INKEY?: IF Q?=" " THEN GOTO 8360 
IF Q?="Y" OR Q?="y" THEN GOTO 8420 
GOTO 8360 
PRINT TIME? 

REM ***************************************** 

REM * GET THE NUMBER AND STORE IN THE * 

REM * FILE * 

REM ***************************************** 

DSEW?="C: C" 

FOR 1=1 TO ND 

IF I<10 THEN DSFILE?=DSEW?+RIGHT?(STR?(I),l)+"-WFM" 

IF I>9 THEN DSFILE?=DSEW?+RIGHT?(STR?( I) ,2)+". WFM" 

REM 

REM * **************************************** 

REM * DIGITIZE CHANNEL 4 * 

REM * 

VEW?=": DIGITIZE CHANNEL4" 

REM 

GOSUB 5840 

VEW?=": WAVEFORM: SOURCE WMEM0RY4; FORMAT ASCII" 

GOSUB 5840 

VEW?=": WAVEFORM: DATA?" 

GOSUB 5840 
REM 

CALL IBRDF(DSO%,DSFILE?) 

CLOSE# 1 
NEXT I 

REM ***************************************** 

REM * FINISH DIGITIZING * 

REM ********* ****** ************************** 

PRINT TIME? 

DSFILE?="HEADER" 

VEW?=": WAV: PRE?" 

GOSUB 5840 
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9180 CALL IBRDF(DSO%,DSFILE$) 

9200 PRINT "HEADER COMPLETE" 

9220 GOSUB 10200 

9240 PR I NT 1 T ********** *********** 1 1 

9260 PRINT "WOULD YOU LIKE TO GET A HARD COPY?(Y/N)" 

9280 Q$=INKEY$: IF Q$=" " THEN 9280 
9300 IF Q$="N" OR Q$="n" THEN GOTO 10080 

9320 IF Q$="Y" OR Q$="y" THEN GOTO 9360 

9340 GOTO 9280 

9360 REM **************************************** 

9380 REM * PLOT THE LAST WAVEFORM * 

9 40 0 REM **************************************** 

9420 VEW$=": HARD: SOUR FACT,WMEM4" 

9440 GOSUB 5840 

9460 PRINT "PLOT ? (Y/NV' 

9480 Q$=INKEY$: IF Q$=" 0 THEN 9480 
9500 IF Q$="N" OR Q$="n" THEN GOTO 9600 

9520 IF Q$="Y" OR Q$="y" THEN GOTO 9560 

9540 GOTO 9480 
9560 VEW$=": PLOT?; " 

9580 GOSUB 5840 
9570 GOTO 9780 

9600 REM ***************************************** 

9620 REM * PRINTING * 

9640 REM ***************************************** 

9660 VEW$=": HARDCOPY: PRINTER DEFAULT" 

9680 GOSUB 5840 

9700 VEW$=": HARDCOPY: PAGE AUTOMATIC" 

9720 GOSUB 5840 
9740 VEW$=": PRINT?; " 

9760 GOSUB 5840 
9780 BOAD$="GPIBO" 

9800 CALL IBFIND(BOAD$,BOA%) 

9820 CALL IBSIC(BOA%) 

9840 VEW$="?Q! " 

9860 CALL IBCMD (BOA%,VEW$) 

9880 V%=1 

9900 CALL IB GTS (BOA%,V%) 

9920 PR I NT 1 1 ************************************ ** 1 1 
9940 PRINT "PLEASE WAIT UNTIL THE PLOTTER HAD FINISHED" 
9960 PRINT "PLEASE RELEASE THE DPO TO LOCAL(Y)" 

9980 Q$=INKEY$: IF Q$=" " THEN 9280 
10000 IF Q$="Y" OR Q$="y" THEN GOTO 10040 
10020 GOTO 9980 
10040 CALL IBLOC(DSO%) 

10060 CALL IBLOC(BOA%) 

10080 RETURN 

10100 REM **************************************** 

10120 REM SUBROUTINE READ DATA TO STRING REF. GURU D-47 
10140 CALL IBRD(DSO%,RD$) 

10160 RETURN 

10180 REM **************************************** 

10200 REM * PROGRAM TO ARRANGE FORMAT FILE * 

10220 REM * SCALE THE DATA * 

10240 REM * * 

10260 REM * * 
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10280 

10300 

10320 

10340 

10360 

10380 

10400 

10420 

10440 

10460 

10480 

10500 

10520 

10540 

10560 

10580 

10600 

10620 

10640 

10660 

10680 

10700 

10720 

10740 

10760 

10780 

10800 

10820 

10840 

10860 

10880 

10885 

10890 

10892 

10896 

10897 

10898 

10899 

10900 

10901 

10902 

10903 

10904 



D% 

E 

F 

G% 

H 

IY 

J 

K 



REM **************************************** 
OPEN "C: HEADER" FOR INPUT AS#1 
INPUT# 1 , A% , B% , C% , D% , E , F , G% , H , I , J , K 
PRINT "FORMAT ";A% 

PRINT "TYPE ";B% 

PRINT "POINT "; C% 

PRINT "COUNT 
PRINT "XINCRE 
PRINT "XOR 
PRINT "XREF 
PRINT "YINC 
PRINT "YOR 
PRINT "YREF 
PRINT "RANGE 
CLOSE# 1 
REM DSE$="C 
FOR K=1 TO ND 
IF K<10 THEN DSFILE$=DSEW$+RIGHT$(STR$(K) , 1)+ 
IF K>9 THEN DSFILE$=DSEW$+RIGHT$(STR$(K) , 2)+ 
IF K>99 THEN DSFILE$=DSEW$+RIGHT$(STR$(K) , 3)+ 
OPEN DSFILE$ FOR INPUT AS #1 
REM 
REM 
REM 
REM 
REM 
REM 

FOR 1=1 TO C% 

INPUT # 1 , WFM 
W(I)“ ( (WFM-J)*H+IY) 

NEXT I 
SUMM=0. 0 

FOR 1=1 TO C% 

SUMM=SUMM+W( I) 

NEXT I 

REM *************************************** 
REM * REMOVE MEAN * 

REM *** 

SMEAN=SUMM/C% 

PRINT "MEAN ="; SMEAN 
FOR 1=1 TO C% 

RRR( I ) =W( I ) - SMEAN 
NEXT I 



wfm" 

wfm" 

wfm" 



10910 CLOSE# 1 

10920 DSFILE$=DSE$+". WFM" 

10940 OPEN "A" ,#1 ,DSFILE$ 

10960 PRINT# 1,T$ 

10980 PRINT# 1 ,W$ 

11000 PRINT# 1,D$ 

11020 PRINT# 1 , DUMMY $ 

11040 PRINT# 1,D% 

11060 PRINT#1, USING "###.##"; FI 
11080 PRINT# 1,ND 
11100 PRINT# 1, PI 
11120 FOR 1=1 TO C % 

11140 PRINT# 1, USING"##.##### "jRRR(I) 
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11160 

11180 

11200 

11220 

11221 

11240 

11260 

11280 

11300 

11320 

11340 

11360 

11380 

11400 

11420 

11440 

11460 

11480 

11500 

11520 

11540 

11560 

11580 

11600 

11620 

11640 

11660 

11680 

11700 

11720 

11740 

11760 

11780 

11800 

11820 

11840 

11860 

11880 

11900 

11920 

11940 

11960 

11980 

12000 

12020 

12040 

12060 

12080 

12100 

12120 

12140 

12160 

12180 

12200 

12220 

12240 



NEXT I 
REN 

CLOSE# 1 
NEXT K 
BEEP: BEEP 
RETURN 



REM 


* 


PI 


= 


NUMBER OF POINTS 


* 


REM 


* 


N 


= 


NUMBER OF SUB AVERAGES 


* 


REM 


* 


ND 


= 


NUMBER OF DATA BLOCKS 


* 


REM 


* 


T$ 


= 


TARGET IDENTIFICATION 


* 


REM 


* 


D$ 


= 


DATE OF MEASUREMENT 


* 


REM 


* 


T2 


= 


DPO SAMPLING INTERVAL 


* 


REM 


* 


W$ 


= 


WAVEFORM TYPE 


* 


REM 


* 


FI 


= 


DPO TIME WINDOW 


* 


REM 


* 


SI 


= 


DPO MAXIMUM SCALE 


* 


REM 


* 


DSE$ 


zz 


OUTPUT DATA FILE NAME 


* 


REM 


* 


DUMMY$ 


= DUMMY $ 


* 



REM 

TS$=SPACE$( 64) 

DUMMY$=SPACE$( 64) 

KEY OFF: SCREEN 0,0,0 

OPEN "DATA. DAT" FOR INPUT AS #1 

INPUT#1 ,P1: INPUT, n,N: INPUT# 1,ND: INPUT#1,T$: INPUT#1,D$: INPUT# 1,T2 
INPUT#1,W$: INPUT# 1, FI: INPUT#1,S1: INPUT#1,DSE$: I NPUT# 1 , DUMMY $ 
INPUT#1 , DELAY 1: CLOSE# 1 
COLOR 7,1,3: CLS 
X=4: Y=16 



PRINT"PROGRAM DATA ACQUISTION " 

LOCATE X,Y+15: PRINT"PRESS THE FOLLOWING KEY" 
LOCATE X+2,Y: PRINT " 

LOCATE X+3,Y: PRINT " 

LOCATE X+4,Y: PRINT " 

LOCATE X+5,Y: PRINT " 

LOCATE X+6,Y: PRINT " 

LOCATE X+7,Y: PRINT " 

LOCATE X+8,Y: PRINT " 

LOCATE X+9,Y: PRINT " 

LOCATE X+ 10, Y: PRINT 
LOCATE X+11,Y: PRINT 
LOCATE X+12 , YfPRINT 
LOCATE X+13,Y: PRINT 
LOCATE X+14,Y: PRINT 
Q$=INKEY$: IF Q$=' 

IF Q$="G" OR Q$= 



F 

P 

A 

B 

T 

D 

W 

K 



< 

< 

< 

< 

< 



S 

R 

U 

X 

G 



FILE NAME FOR DATA OUTPUT. 

NUMBER OF POINTS 

NUMBER OF SUB AVERAGES 

NUMBER OF DATA BLOCKS 

TARGET ID. ($) 

DATE ( $ ) 

WAVEFROM TYPE 

DPO TIME WINDOW 

DPO SAMPLING INTERVAL 

RUN THE PROGRAM 

DPO MAXIMUM VERT MV 

AUTO TIME SET 

DELAY TIME 



DSE$ 

PI 

N 

ND 

T$ 

D$ 

W$ 

FI 

T2 

SI 

DUMMY $ 
38 



IF Q$="F" OR Q$=' 
IF Q$="P" OR Q$=' 
IF Q$="A" OR Q$=' 
IF Q$="B" OR Q$=’ 
IF Q$="T" OR Q$=' 
IF Q$="D" 

IF Q$="W" 

IF Q$="S" OR Q$=’ 
IF Q$="K" OR Q$=' 
IF Q$="R" OR Q$=' 
IF Q$="U" OR Q$=* 



OR Q$=* 
OR Q$= 



g 



THEN 12000 
THEN 12900 



f" THEN 12300 
p" THEN 12340 
a" THEN 12560 
b" THEN 12600 
t" THEN 12640 
d" THEN 12680 
w" THEN 12700 
s" THEN 12740 
k" THEN 12780 
r" THEN 4040 
u" THEN 12820 
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12260 

12280 

12300 

12320 

12340 

12360 

12380 

12400 

12420 

12440 

12460 

12480 

12500 

12520 

12540 

12560 

12580 

12600 

12620 

12640 

12660 

12680 

12700 

12720 

12740 

12760 

12780 

12800 

12820 

12840 

12860 

12880 

12900 

12920 

12940 

12960 

12980 

13000 

13020 

13040 

13060 



IF Q$="X" OR Q$="x" THEN 12860 



GOTO 11720 

PRINT" INPUT FILE NAME( $) EX. .C:F" 
INPUT DSE$: GOTO 12940 



PRINT "A). 



,8: PRINT"NUMBER 


OF POINT ' 


1024 


. . . POINT. ' 


. 512 


. . . POINT. ' 


. 256 


. . . POINT. ' 


. 128 


. . . POINT. ' 



Q$=INKEY$: IF Q$="" THEN 12440 



IF Q$= A OR Q$= a THEN Pl=1024 
IF Q$="B" OR Q$="b" THEN Pl=512 
IF Q$="C" OR Q$="c" THEN Pl=256 



IF Q$="D" OR Q$="d" THEN Pl=128 



GOTO 12940 

PRINT" INPUT NUMBER OF SUB -AVERAGE" 

INPUT N: GOTO 12940 

PRINT "INPUT NUMBER OF DATA BLOCK" 

INPUT ND : GOTO 12940 
PRINT" INPUT TARGET ID($)" 

INPUT T$: GOTO 12940 

PRINT"SET DATE": D$=DATE$: GOTO 12940 

PRINT" INPUT WAVEFROM TYPE($V* 

INPUT W$: GOTO 12940 

PRINT" INPUT DPO SAMPLING INTERVAL" 

INPUT T2: GOTO 12940 
PRINT"DPO TIME WINDOW " 

INPUT FI: GOTO 12940 

PRINT" INPUT DPO VERTICAL SCALE" 

INPUT SI: GOTO 12940 
PRINT" AUTO TIME" 

DUMMY$=TIME$: GOTO 12940 
PRINT" INPUT TIME DELAY" 

INPUT DELAY 1: GOTO 12940 

OPEN "DATA. DAT" FOR OUTPUT AS#1 

PRINT#1 ,P1: PRINT#1 ,N: PRINT#1 ,ND: PRINT#1,T$: PRINT#1,D$ 
PRINT#1,T2: PRINT#1,W$: PRINT#1 , USING "###. FI: PRINT#1,S1 
PRINT#1,DSE$: PRINT#1 ,DUMMY$: PRINT# 1 , DELAY 1 
CLOSE# 1 : CLS 
REM 

GOTO 11700 
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APPENDIX C. CHARACTERISTICS OF LOW LOSS CABLE 



Table 21. CHARACTERISTICS OF PORT RETURN CABLES 



Impedance 


50 ohms 


Capacitance 


26 pF/ft 


Time delay 


1.2 ns/ft 


Velocity of propagation 


85% of light velocity 


Jacket withstand 


1.0 kV 


Center conductor 


15 AWG solid 


Minimum bend radius 


1 inch 


Dielectric constant 


1.4 


Dielectric withstand 


1.0 kV 



Table 22. INSERTION LOSS 



0-4 GHz 


less than 0.36 dB 


4-8 GHz 


less than 0.5 dB 


8-12 GHz 


less than 0.6 dB 


12-16 GHz 


less than 0.72 dB 


16-18 GHz 


less than 0.8 dB 


18-26.5 GHz 


less than 1.0 dB 
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